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Abstract 

Changes in the mean square nuclear charge radii along the lithium isotopic chain were deter- 
mined using a combination of precise isotope shift measurements and theoretical atomic structure 
calculations. Nuclear charge radii of light elements are of high interest due to the appearance of 
the nuclear halo phenomenon in this region of the nuclear chart. During the past years we have 
developed a new laser spectroscopic approach to determine the charge radii of lithium isotopes 
which combines high sensitivity, speed, and accuracy to measure the extremely small field shift of 
an 8 ms lifetime isotope with production rates on the order of only 10 000 atoms/s. The method 
was applied to all bound isotopes of lithium including the two- neutron halo isotope 11 Li at the 
on-line isotope separators at GSI, Darmstadt, Germany and at TRIUMF, Vancouver, Canada. 
We describe the laser spectroscopic method in detail, present updated and improved values from 
theory and experiment, and discuss the results. 
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I. INTRODUCTION 



Laser spectroscopy of lithium isotopes has recently attracted much interest. This is 
for two reasons: First, as a three-electron system it can be used to test the fundamental 
theoretical description of few-electron systems at high accuracy. Second, the extraction of 
nuclear charge radii from isotope shifts for very light elements became possible by combining 
high-accuracy measurements with atomic theory. These charge radii are of special interest 
since one of the lithium isotopes, n Li, is the best investigated halo nucleus, but a nuclear- 
model-independent value of its charge radius was not known until recently. Our ToPLiS 1 



Rial. 



Here we 



Collaboration succeeded in measuring the charge radii of all lithium isotopes 
will describe the experimental setup in detail and discuss the results, including updated and 
improved values from theory and experiment. 

Laser spectroscopy has considerably contributed to our knowledge of ground-state prop- 
erties of short-lived isotopes. Information on nuclear spins, charge radii, magnetic dipole 
moments and spectroscopic electric quadrupole moments can be extracted from isotope shift 
and hyperfme structure measurements in atomic transitions. A particular strength of these 
methods is that theyprovide nuclear-model-independent data. This field has been regularly 
reviewed, see, e.g., 4K3] within the last decades. While nuclear moments can be obtained 
with laser spectroscopy from the lightest to the heaviest elements, it was so far not possible 



to determine nuclear charge radii for short-lived isotopes lighter than neon 

The reason is that the isotope shift, in which the charge radius information is encoded, 
has two contributions: One is the change in nuclear mass between two isotopes (mass shift, 
MS) and the second one is the difference in the charge distribution inside the nucleus (volume 
shift or field shift, FS). The MS by far dominates the isotope shift in light elements and then 
decreases rapidly approximately with increasing mass number A by A~ 2 . By contrast, the 
FS is more than 10 000 times smaller than the MS in the case of the lithium isotopes but 
increases in proportion to Z 2 A~ 1 ^ 3 with the nuclear charge number Z and by far exceeds the 
MS in heavy elements. Separating the tiny FS in the light elements is a complicated task 
and could only be performed on very simple and stable atoms or ions containing not more 
than two electrons. The general approach is a very accurate calculation of parts contributing 



1 Two-Photon Lithium Spectroscopy 
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to the mass shift and the comparison with the experimentally observed isotope shift, as first 
demonstrated in 1993 for the case of helium by Drake [k]]. The difference between these 
values is then attributed to the change in the nuclear charge radius. Consequently, relative 
accuracy better than 10~ 5 must be reached in experiment as well as in atomic structure 
calculations. Here, the correlations in the atom between more than two electrons were an 
unresolved task. Therefore until a few years ago, this was only possible for one- and two 



electron^, and this approach was used for the stable isotopes of hydrogen 



11 



helium 
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13|, and Li+ [14j. 

The extension of such measurements to short-lived isotopes was a challenge. The exci- 
tation and detection scheme used for measuring isotope shifts in the lithium isotope chain 
must provide both high efficiency, because n Li is produced with production rates of only a 
few thousand atoms/s, and high resolution and accuracy for extraction of the tiny field shift 
contribution. Furthermore, it must be a rapid technique since n Li is the shortest-lived iso- 
tope (T1/2 = 8 ms) that has ever been addressed with cw lasers as required for the necessary 
accuracy 2 . To probe the nuclear charge radius, a transition out of the atomic ground state 
of lithium seems to be the natural choice. Resonance spectroscopy in the triplet system of 
the Li + system [jjj] could in principle be performed as well, but the preparation of such a 
metastable state is usually quite inefficient and therefore not appropriate. Spectroscopy on 
trapped ions or atoms has proved to provide the required sensitivity and accuracy for such 
measurements H]; however, the n Li lifetime is much too short to allow sufficient time 
for trapping and cooling. The 2s 2 Si/2 — > 2p 2 Pj transitions in Li were used previously for 
several investigations of 8 > 9 ' n Li in collinear laser spectroscopy by optical pumping and /?- 



NMR detection 



18 



22|. In these cases the lithium ions, provided by the ISOLDE radioactive 
beam facility at CERN, were neutralized in flight in a charge exchange cell and then used as 
a fast atomic beam. However, the Doppler shift of the transition frequency due to the beam 
velocity of the accelerated ions, and the limited knowledge of the acceleration voltage with 
an uncertainty of typically 10 -4 is prohibitive for isotope shift measurements. Recently, such 
measurements became practicable for Be + ions by using a frequency comb and simultane- 



ous spectroscopy in parallel and antiparallel direction 



23 



24j in order to eliminate voltage 



2 Only some radioisotopes of americium, investigated by the group of H. Backe 15|, have even shorter half 
lives but here broad-band pulsed lasers could be applied. 



4 



uncertainties. Those studies are still limited by statistical and systematic uncertainties to 
an accuracy of about 1 MHz, which is not sufficient for lithium isotopes. Hence, in the 
experiment discussed here, the Li + ions are stopped and neutralized before spectroscopy is 
performed on the atomic cloud. To avoid large Doppler broadening, the 2s 2 Si/2 —> 3s 2 Si/2 
two-photon transition is studied, which is to first order Doppler-free. 



After a series of short letters 



25] in which we reported results on the charge radii 



of all isotopes, we describe here in detail the experimental technique as well as the the- 
oretical calculations and discuss all possible sources of systematic uncertainties that were 
further studied. The paper is organized as follows: In Section [III we give a summary of 
the theoretical mass shift calculations that were performed and steadily improved during 
the last years, including relativistic, QED and nuclear structure corrections. Afterwards, 
the experimental setup is described, including radioactive beam production, separation and 
the transformation into a thermal atomic beam. Resonance ionization mass spectroscopy 
is combined with two-photon spectroscopy to reach the required accuracy and efficiency 
simultaneously. The corresponding results and the change in nuclear charge radius for the 
pair of stable isotopes 6 ' 7 Li are presented in Section IIVI Results for the stable isotopes are 
compared with other experimental data obtained by optical spectroscopy of the lithium ion 
or atom to verify the consistency of the data. The paper closes with the presentation of 
the extracted changes in the mean square (ms) nuclear charge radii 5 (r 2 } 6 ' A relative to the 
reference isotope 6 Li. These values can be extracted independently from further information 
in a model-independent way J4]. 

Most interesting for nuclear-structure studies are absolute nuclear charge or nuclear mat- 
ter radii. To determine charge radii from our isotope shift measurements, we need the 
nuclear charge radius of at least one reference isotope obtained by a different technique. A 
nuclear charge radius determination solely based on optical measurements of the absolute 
iransition frequency and atomic theory calculations is thus far only possible for hydrogen 
261 ] . Experimental data are now in principle also available for the lithium isotopes 271 ]: 
however, atomic theory for three-electron atoms is not yet able to achieve absolute transi- 
tion frequencies sufficiently accurate to extract the 10 -9 contribution of the finite nuclear 
size effect. The choice of the reference isotope and the respective charge radii together with 



an evaluation of theoretical results will be discussed in a following publication 



28 
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II. THEORY 



To extract nuclear charge radii from isotope shift measurements of the lightest elements, 
calculations of the mass-dependent part of the isotope shift at the highest level that is 
currently achievable are required. In this Section we will present the latest results for the 
lithium isotopes. The calculation starts with constructing the solution of the nonrelativistic 
Schrodinger equation using variational calculations in Hylleraas coordinates. These solutions 
are then the basis for the calculation of relativistic and quantum electro dynamic contribu- 
tions. Finally, nuclear structure corrections are included. The finite nuclear size effect is 
actually the part that will be extracted from experiment. But the proportionality coefficient 
between the change in the mean square nuclear charge radius and the extracted nuclear 
volume effect in the isotope shift has to be provided by theory. Additionally, as with many 
halo nuclei, the 11 Li nucleus possesses strong electric dipole transitions to low-lying states, 
making it a so-called "soft-dipole" with a relatively large nuclear polarizability compared 
with less exotic nuclei. The influence of this polarizability on the atomic electron levels was 



recently evaluated 



29| and leads to a significant contribution. In the following Sections, the 



details of the atomic structure calculations will be explained. 



A. General Approach 

In order to extract the rms nuclear charge radius from the measured isotope shift, we 
begin by writing the isotope shift for an atomic transition a — > b between isotopes A and A' 
in the form 

r A,A' A' A /i\ 

= 6uSU(A, A') + C a ^ b [fl(A') - fl{A)\ 

r A,A' r A, A' . n r / 2\A,A' /r>\ 

where f c (A) and f c (A') are the rms charge radii for the two isotopes, and the remaining 
term 5v^\ h {A, A') comes from the mass dependence of the atomic energy levels. Hence, it 
is called the mass shift (MS). As will be seen, the coefficient C a ^b is nearly independent of 
the isotopes involved, but in the case of lithium, relativistic corrections to the wave function 



at the origin and leading recoil corrections should be included 30(. The main challenge is 
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to calculate the term 5v a \ h (A,A') to sufficient accuracy. Equation ([I]) can be written in 
the form ([2]) that usually appears in experimental papers, where the particular transition 
index is suppressed. Here, the mass-shift and the 5 (r 2 ) A ' A dependence of the field shift are 
clearly stated. In the following, we will always refer to the 2s 2 Si/2 — > 3s 2 Si/2 two-photon 
transition. 

To give an overview of the contributions to Sv^s , it is convenient to arrange them in 
the form of a double series expansion in powers of a ~ 1/137.036 and the electron reduced- 
mass-to-nuclear-mass ratio \ijM ~ 10 -4 , where /i = mM/(m + M). Table d summarizes the 
various contributions, including the QED corrections and the finite nuclear size term. Since 
all the lower-order terms can now be calculated to very high accuracy, including the QED 
terms of order a 3 Ry, the dominant source of uncertainty comes from the QED corrections 
of order a 4 Ry or higher. For the isotope shift, the QED terms independent of fi/M cancel, 
and so it is only the radiative recoil terms of order a 4 fi/M ~ 10 _12 Ry (~10 kHz) that 
contribute to the uncertainty. Since this is much less than the finite nuclear size correction 
of about 1 MHz, the comparison between theory and experiment clearly provides a means 
to determine the nuclear size. This is the key point to keep in mind when considering the 
theoretical contributions to the isotope shift. 



B. Solution to the Nonrelativistic Schrodinger Equation 

The foundation for the calculation, and the subsequent evaluation of relativistic and QED 
corrections, is to find high-precision solutions to the nonrelativistic Schrodinger equation for 
finite nuclear mass. The past 20 years have seen important advances in developing special- 
ized techniques for doing this in the case of the three-body problem (helium-like systems) 
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31 



321 ] . and more recently the four-body problem 29l. I33l-l36l]. The usual methods of the 



oretical atomic physics, such as the Hartree-Fock approximation or configuration interaction 
methods, are not capable of yielding results of spectroscopic accuracy, and so specialized 
techniques of the type used here are needed. 

For convenience, we begin by rescaling distances according to r — > (m//i)r. The advan- 
tage gained is that the Hamiltonian for a three-electron atomic system can then be written 
in the form 



H = H + \H 



MP 



(3) 



TABLE I: Contributions to the energy in units of the Rydberg constant and their leading orders 
of magnitude in terms of Z, fi/M ~ 10~ 4 , and a ~ 10~ 2 . ad,nuc is the nuclear dipole polarizability 
and ao the Bohr atomic radius. 

Contribution Magnitude 



Nonrelativistic energy Z 2 

Mass polarization Z 2 fi/M 

Second-order mass polarization Z 2 (fi/M) 2 

Relativistic corrections Z 4 a 2 

Relativistic recoil Z 4 a 2 /j,/M 

Anomalous magnetic moment Z 4 a 3 

Hyperfine structure Z gj/J,Q 

Lamb shift Z^o? In a + • • • 

Radiative recoil Z 4 a 3 (ln a)fJ,/M 

Finite nuclear size Z 4 (f c /ao} 2 

Finite size recoil Z 4 [i/M {r c /ao) 2 

Nuclear polarization Z 3 e 2 ; a<i nuc / '(aao) 



with 

3 

H 



1 = 1 1=1 l>] J 



and 

3 

H MP = ^Vi-Vj-, (5) 

in units of 2Rm, where the Rydberg constant for finite nuclear mass is denned by Rm = 
(1 — fj,/M)R 00 , and A = —fi/M can be treated as a perturbation parameter. The Schrodinger 
equation 

HV = Em (6) 
was solved perturbatively by expanding \1> and E according to 

* = *o + A*i + -- - , (7) 

E = e + Xe 1 + X 2 e 2 H . (8) 

8 



Thus Eq. ([6]) becomes 

H V = ^o, (9) 
(£0-^0)^1 = (#mp - £i)*o ■ (10) 

ei and 62 are 

£l = (# |#Mp|*o), (11) 

£ 2 = (* |^Mp|*i)-ei(*o|*i). (12) 

Both \l/ an d ^1 were solved variationally in multiple basis sets in Hylleraas coordinates 
containing terms of the form 

r jl r h „J3 „J12 „j23 „J31 „-apr , i-/3 p r2-7pr 3 
'1 '2 '3 '12 '23 '31 e A 

y\hii)ti2,iz ( r ij r 3) xi > (13) 

where yh^f 2 y 12 1 3 is a vector-coupled product of spherical harmonics for the three electrons 
to form a state of total angular momentum L, and Xi is a spin function with spin angular 
momentum 1/2. 



As described previously 



37 



3s 



the basis set is divided into five sectors with different 
values of the scale factors a p , (3 P , and 7 P in each sector, as labelled by the subscript p. The 
15 independent scale factors are then optimized by a global minimization of the energy. 
Except for restrictions on which correlation terms are included in each sector (see Ref. 38] ), 
all terms from (1T3"]) are included such that 



ji + h + h + 3n + J23 + J31 < 



(14) 



and the convergence of the eigenvalues is studied as Q is progressively increased. Further 
details may be found in Ref. [38]. Since Eq. flSJ) is expressed in units of (1 + A) 2i?oo, the 
explicit mass- dependence of E in units of 2/?^ is 

e + \{e + £l ) + \ 2 {e 1 + e 2 ) + 0(A 3 ). (15) 



E 



C. Relativistic and Relativistic Recoil Corrections 



The lowest-order relativistic corrections of 0(a 2 ) and the spin-dependent anomalous mag- 



netic moment corrections of 0(a 3 ) can be written in the form 



31 



391 ] (in atomic units) 



E rcl = WH^Wj. 



(16) 



where \I/ is a nonrelativistic wave function and H rc \ is defined by 



In (HID, 



H ve \ = B 1 + B 2 + -83c + B,3 Z + -B5 

-7ia 2 ^T (l + -Si ■ sAs^Tij) + 

i>j \ ' 

iw^5(r i ) + ^(A 2 + A 3 ,) 



2 

3"^ ' 3 



-7 ( 2£?3 Z + ^ + -B$ + 2£? 5 



m r 
+7 — A3, . 



(17) 



5! = _:_(V 4 + V 4 + v 4 ^ 



37 ' 



(18) 



5, 



or 



E 



(19) 



2 3 1 

„ at \ 1 , . , 

4 = 7^ T r i* x Pi ' ( s * + 2s i) 

1 i±i ^ 



(20) 



Za 2 ^ 1 

2 tr r ' 



(21) 



3 r 



B 5 = a*J2 



i>j 



r 3 \ S i ' S j> 

ij 



(22) 



r iZa 2 v-^ 

^ E 

j"=i 



1 „ 1 , 



(23) 



3 1 

A 3Z = Za 2 ^ — T i x P ■ s i 



i=l 



(24) 



n 3e 



2 3 1 
9 ^ 3 F ^ X Pi ' S -?' y 



(25) 
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with p = pi + P2 + P3, and 7 is 

7 = ^ + (-0.32847) (-Y + ---. (26) 

The operator —ira 2 X^i>,(l + | s i ' s j)^( r ij) can replaced by na 2 ^i>j $( r ij) an d the ex- 
pectation value of the spin-spin term B 5 vanishes. The terms proportional to m/M are 
the nuclear relativistic recoil corrections and the terms proportional to 7 are the anomalous 
magnetic moment corrections. 

The perturbing effect of mass polarization on the expectation values of Breit operators 
can be obtained using 

* = * + A (*i-(*i|* ) *()) + ••• , (27) 
where the extra term — (\l/ 1 | 1 I'o} x I / o is added to so that the first two terms of the right 



hand side are orthogonal to each other 40J]. Thus, for a Breit operator A, one has 

(V\A\V) = a + Xa 1 + --- , (28) 

where 

a = (* |^o>, (29) 

and 

ai = 2(* |^|* 1 >-2(*o|*i)(*o|^|*o). (30) 

Furthermore, due to the use of ^-scaled atomic units in Eq. ([3]), the units of in 
Eq. 028p are {^/m) n 2R OQ , where — n is the degree of homogeneity of operator A in the 
three-electron coordinate space such that 

A(/3n, 0r 2 , (3r 3 ) = /TMfa, r 2 , r 3 ) . (31) 

Using 

(£f = (l + Ar«l + nA, (32) 
one has the explicit mass- dependent formula 

(*|A|*> = a + A (no + ai) + 0(X 2 ) , (33) 

in units of 2-Roo. 
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D. QED Corrections 



Until recently, the QED contributions of lowest order a 3 Ry were the major source of 
uncertainty in calculations of atomic energy levels and the isotope shift for atoms more 
complicated than hydrogen. However, tremendous progress has been made in recent years. 
Complete results to lowest order a 3 Ry are now readily obtainable, and higher order correc- 
tions can be estimated in a screened hydrogenic approximation. For a many-electron atom, 
it is convenient to express the total QED shift in the form 



E 



QED 



Ela + Em,i + Era + E Li 



2 • 



(34) 



where E^i is the mass- independent part of the electron-nucleus Lamb shift (the Kabir- 
Salpeter term 4l|), Ema contains mass scaling and mass polarization corrections, 
contains recoil corrections (including radiative recoil), and £l,2 is the electron-electron term 



rese in turn. 



44J , except that 



originally obtained by Araki [42| and Sucher |43|. We now discuss each of t 

The term closely resembles the corresponding hydrogenic Lamb shift 4 
an overall multiplying factor of (S(r)) = Z 3 5l,o/ '(7rn 3 ) for the hydrogenic case is replaced by 
the correct expectation value (X]jLi^( r j)) f° r the niulti-electron case, summed over the N 
electrons. The residual state dependence due to other terms such as the Bethe logarithm 
discussed below is then relatively weak. 

Following the notation of Eq. for the mass polarization corrections, the main electron- 
nucleus part for infinite nuclear mass is (in atomic units throughout) 



E 



LA 



4Z« 3 (5(rO) 



(o) 



ln(Za)" 2 - f3{n 2 L) 



19 
30 



(37raZ)0.765 405 577 



a 



+ - [0.404 17 - (3aZ/4)21.556 85] 

7T 



(Zaf 



4 



ln 2 (Za)" 2 + C 61 (ls x nL) ln(Za) 



+ C m (ls x nL)]} , 
the mass scaling and mass polarization corrections are 



(35) 



E 



MA 



E 



LA 



M(S{Ti))(°) 

+ 3M ^~ A(3[ 



MP 



n 2 L 



(36) 
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and the recoil corrections (including radiative recoil) are given by 

4Z 2 a 3 /i (c^) > (0) 



E 



K.l 



3M 



iln(Z«r 2 -2/^ 2 L)-l 



7 3 
--a(ls x nL) - -(vra)l. 36449 

3 1 
+ -7rZaD 50 (n 2 L) + -a 2 Zln 2 (Za)" 2 

4r Zi 



(37) 



These equations involve contributions to the hydrogenic Lamb shift obtained by many au- 
thors, as summarized by Eides et al. 44j |. The quantity f3(ls x nL) = ln(/c / Z 2 Roo) is the 
two- or three-electron Bethe logarithm, and the two terms 1 — A/3mp(1s x nL) in Eq. (|36|) 
account for the mass scaling and mass polarization corrections to (3(ls x nL) respectively. 
These terms are further discussed below. The term a(ls x nL) is a well-known part of the 



hydrogenic Lamb shift. Its many-electron generalization is given by 



34 



45 



46| 



a(ls x nL) 



Ef=i<^-)> 



-3, 



where 



= ^lim^(r7 3 ( e ) +4vr( 7 + ln6)5(r,)) . 



(38) 



(39) 



7 is Euler's constant and e is the radius of a sphere about rj = that is excluded from the 
integration. 

The orders of magnitude for the other state-dependent coefficients Cqi(1s x nL), 
Cqo(1s x nL), and D 50 (ls x nL) are all estimated from the generic formula 

xX(ls) + X(nL)/n 3 



X(ls x nL) 



(40) 



where X(nL) is the corresponding one-electron coefficient, evaluated directly with the full 
nuclear charge Z for L = and with a fully screened nuclear charge Z — x for the outer 
nL electron for L > [471 ] . The numerical values can be immediately calculated from the 
hydrogenic values, as discussed by Drake and Martin ^2] for helium, and Yan and Drake jsjj 
for lithium. Their contribution to the transition energy is taken to be the QED uncertainty. 

The electron-electron QED shift i?L,2 can similarly be separated into mass-independent 
and mass- dependent parts according to 



E 



L,2 



E 



(0) 
L,2 



m El > 2 + 



(41) 
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where 



<) = a 3 (y ln« + ^) $>(r,)>^ - ^a 3 g<°>, (42) 

^ ' i>j 

and the mass scaling and mass polarization corrections are 

<> = + (ylB<« + ^)E« r «)> (1) 

^ ' i>3 

-ya 3 (gW + ^(5(r.,))(°A . (43) 
Following our notation, the Q^-term for infinite mass is given by 

Q i0) = ^SEW £ ) +47T( T + In (44) 

The term is the correction due to the mass polarization correction to the wave function. 
As a word of explanation, the infinitesimal limiting quantity e has dimensions of distance, 
and so it generates an additional finite mass correction when distances are rescaled for the 
finite mass case according to e — > (/x/m)e. 

The principal computational challenge is the calculation of the Bethe logarithm f3(ls 2 nL) 
in E^i, originating from the emission and re-absorption of a virtual photon, and the finite 
mass correction A/?mp(1s 2 nL) due to mass polarization. The Bethe logarithm is the loga- 
rithmic remainder after mass renormalization and is defined by 

/3(ls 2 nL) = ^= (45) 
|(*o|pi + P2\i)\ 2 (Ei - E ) In \Ei - E \ 



J2\(*o\pi + P2\i)\ 2 (Ei-E t 



0) 



The accurate calculation of (3 (Is 2 nL) has been a long-standing problem in atomic physics. 
This has been solved by use of a discrete variational representation of the continuum in terms 



of pseudostates |48j . The key idea is to define a variational basis set containing a huge range 
of distance scales through multiple sets of exponential scale factors a and (3 that themselves 
span many orders of magnitude. The Bethe logarithm comes almost entirely from virtual 
excitations of the inner Is electron to p-states lying high in the photoionization continuum, 
and so the basis set must be extended to very short distances for this particle. The outer 
electrons are to a good approximation just spectators to these virtual excitations. 

14 



TABLE II: Comparison of Bethe logarithms for lithium and its ions. 



Atom Li(ls 2 2s) Li(ls 2 3s) Li+(ls 2 ) Li+ + (ls) 



ln{k /Z 2 Ry) 2.98106(1) 2.982 36(6) 2.982 624 2.984128 



Table [TT] compares the Bethe logarithms for the two lowest s-states of lithium with those 
for the Li-like ions Li + (ls 2 X S) and Li ++ (ls 2 S). The comparison emphasizes that the 
Bethe logarithm is determined almost entirely by the hydrogenic value for the Is electron 
and is rather independent of the state of excitation of the outer electrons or the degree 
of ionization. In order to make the connection with the hydrogenic Bethe logarithm more 
obvious, the quantity tabulated is ln(ko/Z 2 Ry). The effect of dividing by a factor of Z 2 is 
to reduce all the Bethe logarithms to approximately the same number /3(ls) = 2.984 128 556 
for the ground state of hydrogen. 



E. Nuclear Polarizability 

The interaction of the nucleus with an electromagnetic field can be described by the 
Hamiltonian 

H mt = q A - d • E - M ■ B - § (r 2 ) V • E, (46) 

6 

which is valid as long as the characteristic momentum of the electromagnetic field is smaller 
than the inverse of the nuclear size. Otherwise, one has to use a complete description in 
terms of form factors and structure functions. Under this assumption, the dominant term 
for the nuclear excitation is the electric dipole interaction. This is the main approximation 
of this approach, which may not always be valid. It was checked however that higher order 



polarizabilities are quite small (below 1 kHz) for deuterium 49 



50[ , and this is similar for He, 



Li and Be isotopes. Within this low electromagnetic momentum approximation, the nuclear 



polarizability correction to the energy is given by [29| (in units h = c=l,e 2 = 4ira 



E pol = -ma' ($^ 3 (r a )} (m 3 a pol ), (47) 

a 

where m is the electron mass and the expectation value of the Dirac 5 is taken with the 
electron wave function in atomic units. For hydrogenic systems it is equal to Z 3 /n. In the 

15 



equation above, a po \ is a weighted electric polarizability of the nucleus, which is given by the 
double integral 



16 a 



a 



pol 



r dE^Mm? r — 

Je t e 2 J w 



E 



E 2 + w 2 



(«+!)(«*+!) \K+1 



+ 



(48) 



where k — y\ + 2im/w and E? is the excitation energy for the nuclear breakup threshold. 

The kets \4>n) and \E) denote the ground state of the nucleus and a dipole excited state 
with excitation energy E, respectively. The square of the dipole moment is related to the 
so called B(E1) function by the relation 
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(49) 



in units e 2 fm 2 MeV -1 , which explains the presence of e 2 in the denominator in Eq. (148)) . 

If Et is much larger than the electron mass m, one can perform a small electron mass 
expansion and obtain a simplified formula 5l|, |52[: 
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with the static electric dipole polarizability 
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and the logarithmically modified polarizability 
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This approximation is for example valid for 3 He and 4 He isotopes, but not for n Li or n Be. 

In the opposite situation, i.e., when m is much larger than Et, we can use the nonrela- 
tivistic limit, and the polarizability correction adopts the form (with m being the reduced 
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FIG. 1: (Color online) Electric dipole line strength by Nakamura et al. [17|] adapted to the new 



value of Et = 369.15(65) keV from 
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This approximation is justified for muonic helium atoms or ions, because the muon mass 
(~106 MeV) is much larger than the threshold energy Et = 369.15(65) keV 53|. This 
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formula requires, however, a few significant corrections, namely Coulomb distortion and form 



factor corrections. They were obtained by Friar in |54J for the calculation of the polarizability 
correction in /i 4 He. This nonrelativistic approximation, however, is not valid for electronic 
atoms since the typical nuclear excitation energy in light nuclei is larger than the electron 
mass. This indicates a possible limitation in the comparison of charge radii obtained from 
electronic and muonic atoms, as the nuclear structure effects are much different. 

The B(E1) function is directly related to the photoabsorption cross Section at photon 
energies E 

16tt 3 „/l dB(El)\ 

which allows us to obtain the B(E1) function from experimental data presented in 17J and 
shown in Fig. [TJ With the two-neutron separation energy = 369.15(65) keV from j^j 
one obtains 

a po i = 60.9(6.1) fnT 3 = 1.06(0.11) 1(T 6 ^, (55) 

and a polarizability correction to the 2s 2 Si/ 2 — > 3s 2 Si/ 2 transition frequency of v vo \ = 
39(4) kHz. The polarizability correction for the lighter lithium isotopes is expected to be 
negligible due to much larger separation energies. 



F. Extraction of Nuclear Charge Radii 

The various contributions to the isotope shift as they have been discussed so far are 
listed in Table HTT1 for all lithium isotopes relative to 6 Li. The mass values that were used for 
the calculations are also included for reference. In addition, there is a significant electronic 
binding energy correction AM = — Sbmdmg/c 2 - The terms are classified according to their 
dependence on fi/M and the fine structure constant a. As an important check, most of 
;he results have been calculated independently by Puchalski and Pachucki (P&P) in Poland 



571 ]. and by Yan and Drake (Y&D) in Canada 58|. Two exceptions are the Bethe logarithm 



part of the radiative recoil correction, which have been calculated only by Y&D 



5s 



and 



the nuclear polarizability correction, which has been calculated only by P&P (29(. In some 
cases, the two sets of results are slightly different due to different methods of calculation, 
as noted in the table. In these cases, both sets of results are given for comparison. The 
differences are not large enough to affect the determination of the nuclear charge radius from 
experimental results, but they indicate the areas where further work is desirable. 
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TABLE III: Contributions to the mass shift <^^ms Theory °^ (A = 7, 8, 9, 11) relative to 6 Li 
in the 2s 2 Si/2 —> 3s 2 Si/2 transition. The mass-dependent terms are calculated using the masses 

55j. The 



M listed in the first row. The mass of the reference isotope 6 Li is 6.015122794(16) amu 
nuclear polarizability contribution v po \ is included for n Li and expected to be negligible for the 
other isotopes. The unit of the electronic factor Cq : a is MHz/fm 2 . All other values are in MHz. 
Both sets of theoretical results are given in cases where they differ (see text). 



Term 


7 Li 


8 Li 


9 Li 


n Li 


M (amu) 


7.0160034256(45) a 


8.02248624(12) b 


9.02679020(21) b 


11.04372361(69) b 


H/M 


11 454.655 2(2) c 


20 090.837 3(9)° 


26 788.479 2(13)° 


36 559.175 4(27)° 


(m/m) 2 


-1.7940 


-2.9644 


-3.764 2 


-4.761 9 


a 2 fi/M 


0.017 2 d 


0.030 2 d 


0.040 2 d 


0.055 d 




0.016 8(1)° 


0.029 5(2)° 


0.039 3(3)° 


0.053 7(4)° 


a 3 fi/M 


-0.048 5(6) 


-0.0851(11) 


-0.113 5(15) 


-0.154 8(21) 


a 4 fi/M 


-0.009 2(23) d 


-0.016 l(40) d 


-0.021 5(63) d 


-0.029 4(73) d 




-0.008 4(28)° 


-0.014 7(41)° 


-0.019 6(66) e 


-0.026 8(90)° 


^pol 








0.039(4) 


Total 


11 452.820 7(24) d 


20 087.801 9(42) d 


26 784.620 2(66) d 


36 554.323(9) d 




11452.821 1(28)° 


20 087.802 6(50)° 


26 784.6213(67)° 


36 554.325(9)° 




-1.571 9(16) f 


-1.571 9(16) f 


-1.572 0(16) f 


-1.570 3(16) f 



Ref. 



53| 



° Uncertainties for this line are dominated by the nuclear mass uncertainty. 
d Calculation by Puchalski and Pachucki (this work). 
e Calculation by Yan and Drake (this work). 

A 25% error is assumed for the relativistic correction due to the estimation of the relativistic 
correction to the wave function at the origin on the basis of a known result for hydrogenic 
systems. 
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The various contributions are as follows: The term labelled fi/M contains the sum 
of the reduced mass scaling of the nonrelativistic ionization energy, and the first-order 
mass polarization correction. For example, the mass-scaling term is 2R 00 [E^(Li + , oo) — 
E NR (Li, oo)][(h/M)a u — (/i/M)e Li ], where -Enr^LI, oo) is the nonrelativistic energy for infi- 
nite nuclear mass (in atomic units), and i?oo is the corresponding Rydberg constant. 

The term of order (n/M) 2 comes from second-order mass polarization. The relativistic 
recoil terms of order a 2 fi/M results from mass scaling, mass polarization, and the Stone 
terms as expressed by Eqs. (|23|) and (|24|) . The two sets of results here do not quite agree 
because the Y&D results include partial contributions from the next-higher-order terms 
of order a 2 (fi/M) 2 . The difference between the two calculations is an indication of the 
uncertainty due to the partial neglection of these higher-order terms. The numbers in this 
row are anomalously small because of nearly complete numerical cancelation between the 
mass polarization and mass scaling plus Stone contributions. For example, for the case of 
Li, the two parts are -17.9198(4) MHz and +17.9735 MHz, respectively, resulting in a final 
recoil term of only 0.0537(4) MHz. Because of this cancelation, the percentage uncertainty 
is relatively large. 

The radiative recoil terms of order a 3 fi/M similarly come from a combination of mass 



scaling, mass polarization, and higher-order recoil corrections, as discussed in 29i.l33l436l.l58. 



591 ] . The most difficult part of the calculation is the Bethe logarithm contribution obtained 



in Ref. 



The sum of all these terms as given in the penultimate row in Table II III is the mass- 
dependent part of the isotope shift ^^ms Theory between 6 Li and A Li. It includes the nuclear 
polarizability correction and the adjusted relativistic recoil term as calculated for the first 



time in 29]. The difference between this value and the measured isotope shift must arise 
from the finite nuclear size effect. To extract nuclear charge radii information from this 
value, the electronic factor Ca,a> as defined in Eq. (|2]) is required. 

The contribution of the finite nuclear size effect to the total transition energy can be 
written as 

AE QUC = C ■ f 2 c , (56) 
with the mean square nuclear charge radius f 2 = (r 2 ). The electronic factor C can be 
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TABLE IV: Finite size corrections C i - m,n ^ of the constant C in MHz/fm 2 for the 2s —> 3s transition 
in n Li and 6 Li. 



Coefficient 


n Li 


6 Li 


C(4,0) 


-1.566501 


-1.566501 


C( 6 '°) 


-0.006675 


-0.006641 




0.000231 


0.000424 


C 


-1.572945 


-1.572718 



expanded into a power series of a and (fi/M) 

_ R na m 2 c 3 a 2 



^(4,0) + A (7(4,1) + ^(6,0) + _ ^ ) (57) 

where the C^ m ' n ^ are dimensionless coefficients that refer to corrections on the order of 
ot m (ii/M) n . The leading-order term is proportional to the nonrelativistic wave function at 
the origin 

c m = ^l£ A ^(o)^ (58) 



and C^ 6 ' ^ comes from the relativistic correction to the wave function at the origin 60J 



£(6,0) w -Z 2 \n(Zamf c )C^°l (59) 

Contrary to these terms, C^ 4,1 ) is mass- dependent. It consists of a term originating from 
the mass scaling and one from the mass polarization j30|. For example, the contributions of 
these terms to the C-coefficients of the transition frequency for 6 ' n Li are listed in Table HVl 
As compared with the leading term C^'°\ the C^ 6 ' -' and C^ 4,1 -* corrections are of the order 
1CT 3 to 10~ 4 . Thus, they are significant in the systematic studies of the isotope shift at the 
current level of theoretical uncertainties in (^Ms/Theory- 

Since the correction term C^ 6 ' ^ depends on the nuclear size given by f c , we have to use 



29). 



initial approximate values vao for the charge radii, e.g., those which are obtained in 
The influence of this approximation is negligible. With these values we can obtain the Ca,a' 
coefficient for the isotopes A and A' as 

C A,A' = J! ~2 • 



r A'0 r A0 
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The obtained numerical results for constants Ca,a' for the relevant isotope shifts using this 
relation are included in Table IIII1 

Concerning the accuracy of the calculated mass shift contributions, the uncertainty is 
largely dominated by the nuclear mass uncertainty for 11 Li, as listed in the table. An 
exception is the relativistic recoil term of order a 2 fi/M. Because of the almost complete nu- 
merical cancelation of individual contributions that were already mentioned, the percentage 
uncertainty is correspondingly large. 

For convenience, the theoretical result relating the measured isotope shift 5v to the change 
in the mean square nuclear charge radii between 6 Li and n Li is 

^tsXzs = 36 554.324(9) 

+ 1.570 [f 2 ( 6 Li) - f 2 ( u Li)] , (61) 

where frequencies are given in MHz and radii in fm. The results presented in Table II III are 
the basis for the extraction of the change in the mean square charge radius for all lithium 
isotopes by comparison with the experimental values, according to 

[rK A Li) - f 2 ( 6 Li)] = ^-p-^MS,T h eo ry _ (62) 

III. EXPERIMENTAL SETUP 

Figure [2] shows an overall schematic of the experimental apparatus used for laser spec- 
troscopy of the 2s 2 Si/2 —> 3s 2 Si/2 two-photon transition in lithium. Ions of both stable 
and short-lived lithium isotopes, produced at an accelerator and accelerated to 40 keV, are 
mass separated and transported to the experimental area. Here they are stopped and neu- 
tralized in a thin carbon foil catcher. This foil is heated with a CO2 laser beam to about 
1700-1800 °C , so that the implanted lithium atoms diffuse quickly to the surface. Atoms 
released in the forward direction drift into the ionization region in front of a quadrupole 
mass filter (QMF). Here, the lithium atoms are resonantly ionized with laser light at 735 
nm and 610 nm according to the three-step four-photon resonance ionization scheme 

2o 2x735 nm 2 o decay 2 j~, r r r>\ 

2s S1/2 — > — > 3s S1/2 > 2p P 3 / 2 ,i/2 (63) 

rfs30ns 

610 nm „ , 2t-> 610 nm,735 nm T . i 

> 3d D 3/2i 5/2 > Li , 
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with a two-photon transition followed by spontaneous decay of the 3s state with a lifetime of 
t w 30 ns and subsequent resonance ionization as discussed in more detail in Section IIII CI 
For brevity, if the meaning is clear, the laser-driven resonance transitions are abbreviated 
as the 2s — > 3s and 2p — > 3d transitions, respectively. The photo-produced ions are then 
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FIG. 2: (Color online) Experimental setup to measure the 2s 2 Si/2 —> 3s 2 Si/2 electronic transition 
in lithium. CDEM: Continuous dynode electron multiplier, QMF: Quadrupole mass filter, PZT: 
Piezoelectric transducer. 
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mass analyzed with the QMF and detected with a continuous dynode electron multiplier 
(CDEM) detector. The isotope shift in the 2s — > 3s transition is measured by tuning the 
735 nm light across the lithium two-photon resonances. The individual parts of the system 
will be described in detail in the following Sections. 



A. Production of Radioactive Lithium Isotopes 

Radioactive lithium isotopes were produced at the on-line mass separator 3 at GSI Darm- 
stadt in 2003 [l| and in a second experiment at the ISAC mass separator facility at TRIUMF 
in September and October 2004 (see Table |V|) . The experimental setups were almost iden- 
tical for the two experimental runs. If changes were made, the version used for the TRIUMF 
experiment is presented in this publication. 

At GSI, 8 ' 9 Li were produced by directing a 12 C beam from the UNILAC with an energy of 
11.4 MeV/u onto a 100 mg/cm 2 tungsten target. Fast reaction products entering the hot ion 
source of the mass separator through a tungsten window were stopped in a sintered graphite 
catcher. Atoms diffusing out of the catcher were surface- ionized and extracted through a 
small hole. Ion yields towards the experiment are listed in Table |V] 

At TRIUMF, 8 > 9 ' 11 Li were produced with a 500 MeV primary proton beam of 40 \xA ex- 
tracted from the H~ cyclotron. A stack of tantalum foils was used as target in order to allow 
fast release of the short-lived lithium isotopes that are produced by target fragmentation. 
The reaction products were surface ionized and extracted from the ion source with a beam 
energy of 40 keV. Mass sep aration was obtained with a 60° preseparator magnet followed by 
a 120° main separator 6l|. A fast switch (kicker) installed behind the main mass separator 
was used to turn the ion beam on and off. 

The mass separated ion beam was transported into the ISAC low energy experimental 
area, where the ToPLiS experiment was installed. The existing low energy beamline was 
extended to allow installation of deflector plates and a quadrupole doublet, as depicted in 
Fig. [3j for precise shaping and steering of the ion beam. A similar ion optics was existing at 
the GSI mass separator. To detect the small number of n Li ions delivered to the experiment, 
a channeltron-type detector was installed on a linear feedthrough that also carries the carbon 



The mass separator was shutdown in early 2004. 
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FIG. 3: (Color online) Schematics of the ion optical setup (left) in front of the catcher holder 
which holds a channeltron- type continuous dynode multiplier detector for aligning the n Li ion 
beam and several carbon foils of different thickness. Details of the catcher holder are shown on 
the right. Dimensions are given in mm. The holder, fixed to a linear feedthrough, can be moved 
in the x-direction horizontally and vertically in the y-direction. 



TABLE V: Half-lifes and typical ion yields in ions/s for the radioactive lithium ions at GSI and 
ISAC-TRIUMF. 



Beamtime 



^Li 



3 Li 



ii 



Li 



Half-life (ms) 
Yield (GSI 2003) 
Yield (TRIUMF 09/2004) 
Yield (TRIUMF 10/2004) 



838(6) 

3.6xl0 5 

10 8 

8xl0 8 



178.3(4) 
1.8xl0 5 
10 6 

9xl0 7 



8.59(14) 



3.5xl0 4 



catcher foils and could be moved into the beam in front of the QMF. 



B. Neutralization and Atomic Beam Generation 



After production and mass separation, the 40 keV ion beam had to be converted into a 
thermal beam of neutral atoms. This process must be efficient and considerably faster than 
the 8.6 ms half-life of the n Li ions. Therefore, the ion beam was directed onto a thin carbon 
foil (Fig. [3]). In this 'catcher' foil the ions were stopped and neutralized. "Stopping and 
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Range of Ions in Matter" (SRIM) calculations 62J were used to estimate the thickness of foil 
that would stop the ions shortly before they reached the back side of the foil. This enabled 
rapid release of the neutralized atoms in the preferred direction towards the ionization region 
(Fig. [2]). Catcher foils with thickness around 300 nm were prepared in the GSI Target 
Laboratory and glued onto a stainless steal holder, providing three positions equipped with 
catcher foils of different thickness (ranging from about 60 to 80 /ig/cm 2 ). The target holder 
was mounted on a linear feedthrough that moves the holder along the x-position so that 
foils of different thickness could be located in front of the quadrupole mass filter. The holder 
was also movable in y and z directions to optimize the vertical position of the catcher foil 
and the distance to the QMF ion optics. As mentioned before, fast release is essential for 
measurements on 11 Li. In order to achieve this, a 2 mm diameter spot on the catcher foil 
was heated to about 1800°C with a CO2 laser beam that was coupled into the system along 
the axis of the QMF (see below). As shown in Fig. HJ a release of the implanted atoms 
within 0.5 ms was observed under these conditions. The temperature of the catcher foil was 
still sufficiently low that only a very small fraction of about 10~ 4 of the incoming ions were 
released as surface ions. The release time of the lithium atoms from the graphite catcher foil 
was measured by observing the increase of surface ions as a function of time after turning 
the ion beam on. For this purpose the fast kicker at the mass separator was used, the 
electrostatic ion optics of the QMF were set for surface ion detection and the output of 
the CDEM detector was fed into a fast timing amplifier for linear pulse amplification. The 
amplified signal was recorded with a multichannel scaler. The measurement was performed 
with a beam of stable 6 Li and the observed signal with a time resolution of 1 ms and 0.05 
ms is shown in Fig. H](a) and|4](b), respectively. The increase in ion beam intensity from a 
practically background-free baseline is very fast and is not resolved with the time resolution 
of 1 ms. The graph with the higher time resolution (Fig. H](b)) shows two contributions: 
The ion signal rises to about 70% of the final level within 100 /xs and is then followed by 
a slightly slower increase to the full intensity which is reached after about 300 /is. A large 
contribution of ions that have sufficient energy for simply penetrating the foil was never 
observed. Hence, it is assumed that even the flank was originating from ions that were first 
stopped in the catcher foil and released again. However, even if the 300 /xs time scale would 
be the relevant one for the release time of neutral atoms, this amounts to only 3% of the 
n Li half-life and is clearly sufficient to have a good release efficiency. The total conversion 
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FIG. 4: (Color online) Time-resolved surface-ion release recorded for 6 Li to determine the release 
time at a catcher temperature of about 1700-1800° C. The t = position is approximately the time 
when the ion beam was turned on. Time resolution was 0.5 ms in (a) but after the first three data 
points only every tenth channel is shown as an individual data point for clarity. In (b) a resolution 
of 50 /is was used. 

efficiency of the foils was estimated to be about 50 % and the transport efficiency into the 
laser beam to 20%. 



C. Excitation and Ionization Scheme 



The laser excitation and ionization scheme has to provide both, high resolution and high 
efficiency in order to detect the signal for ion yields of only a few thousand ions/s with an 
accuracy of about 0.1 MHz. 

Figure [5] shows the level scheme for neutral lithium with the excitation path chosen for 
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FIG. 5: (Color online) Resonance ionization scheme for lithium as used in the reported experiments. 

resonance ionization: Lithium atoms in the 2s 2 Si/2 ground level are excited via the Doppler- 
free two-photon transition at 735 nm to the 3s 2 Si/2 level. This transition offers a narrow 
resonance and leads to efficient excitation since all velocity classes can be excited simulta- 
neously. However, relatively large intensities are required for saturation. The two-photon 
excitation is followed by a spontaneous decay into the 2p 2 Pj levels. Atoms in the 2p 2 Pj 
levels are subsequently excited into the 3d 2 Dj/ levels with resonant laser light at 610 nm and 
then photoionized by absorption of a photon at either 735 nm or 610 nm. In this way the 
states used for measuring the isotope shift are decoupled from the states used for ionization 
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and detection which results in a strong reduction of the ac Stark broadening and ac Stark 



shifts 



63]. Compared with fluorescence detection, this multi-step excitation and ionization 



followed by mass- selective ion detection has the advantage that the 2s — > 3s resonance can 
be detected with a very high efficiency and an extremely high signal-to-background ratio. 



D. Mass Spectrometer 

To achieve the additional background suppression, mass selective ion detection is used. 
Quadrupole mass filters (QMF) provide ideal conditions to mass separate the ions with 



65] 



thermal energies that are produced in the laser ionization process described above 

Therefore, the laser beams have to be overlapped with the thermal atom cloud in a region 
from which the created ions can be extracted into the QMF. Moreover, surface ions that are 
produced on the hot catcher foil must be efficiently suppressed, because every surface ion 
of the isotope under investigation entering the laser beam cannot be distinguished from a 
photo ion and will produce background events. 

The QMF used for the reported experiment was a commercial instrument from ABB 
Extrel (ABB Extrel, Pittsburgh, USA, Model No. 150 QC) with 9.39 mm radius rods of 
21 cm length and a free-field radius of r = 8.33 mm. The system is driven at a frequency 
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FIG. 6: (Color online) Schematics of the optical resonator, the ion optics, the quadrupole mass 
filter (QMF) and the continuous dynode electron multiplier (CDEM) detector. 
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of 2.9 MHz and can be used for ions up to A/q < 64. This model provides transmission 
close to 100 % and an excellent neighboring mass suppression as has been demonstrated in 
simulations 
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66|, ultratrace analysis applications [67, 
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and measurements discussed 



below. Ions transmitted through the QMF were detected with an off-axis continuous dynode 
electron multiplier (CDEM) detector. The entrance opening of the channeltron was biased 
with about -2000 V such that the positive ions were accelerated in the CDEM. 

For the measurements reported here, the ionization region of the QMF was modified as 
shown in Fig.|6j The original axial electron impact ion source of the EXTREL device was 
removed and the electrodes before the quadrupole structure were replaced by specifically de- 
signed electrodes allowing free access for the laser beams. The distance between the catcher 
foil and the ionization region was 2-3 mm to have an efficient transfer of the released atoms 
into the ionization region. The QMF ion optics could be operated to accept either laser- 
created ions from inside the ion region or surface ions created on the hot carbon catcher 
foil. In the first mode, the ion region (Fig. |6]) was held at 3.9 V relative to the grounded 
catcher foil. This repels catcher surface ions, while neutral atoms can enter the laser ioniza- 
tion region. The potential difference between the catcher foil and the ionization region was 
small enough to ensure that electrons emitted from the hot catcher surface and accelerated 
into the ion region, do not gain sufficient energy to ionize lithium atoms by electron impact 
ionization (Li ionization potential: 5.39 eV). Photo ions were then extracted by the negative 
extractor voltage and focussed into the QMF rod structure using the remaining electrostatic 
lenses. In the second mode, surface ions produced at the hot catcher foil were accelerated 
into the QMF with a negative voltage at the ion region and the extractor was operated as 
another lens for adapting the beam properties to the QMS acceptance. Ion optical settings 
of all lenses for both detection modes are summarized in Table lVll Ions transmitted through 
the rod system were focussed with an exit lens and detected with the CDEM detector model 
DeTech 5402AH-021, which was chosen for its very low dark count rate of typically 5 to 
20 mHz. 

The lithium isotopes were produced with rates that differ by orders of magnitude (see 
Table |V|) . For additional suppression of neighboring isotopes, the QMF was tested with ions 
on the stable isotopes 6,7 Li produced by surface ionization on the hot carbon catcher foil. 
Ion optical parameters of the QMF extraction and transport optics were optimized in order 
to obtain a flat top profile with steep flanks and maximum transmission. An example for the 
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TABLE VI: Settings for the extractor, the electrostatic lenses and the quadrupole mass filter (see 
Fig. [6]) for the detection of either surface ions or photo ions produced by resonance ionization. All 
values are in V. Pole bias is an additional common DC voltage that is applied to all four QMF 
rods to change the kinetic energy of the ions entering the rod system. 
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FIG. 7: (Color online) Mass spectra of the stable lithium isotopes 6 Li and 7 Li. 

peak profile of surface ions is shown in Fig. [7J Typically, a suppression of neighboring masses 
of > 10 8 was achieved. The dark count rate of the CDEM detector was about 10 mHz. 

The entire QMF system consisting of ion source, quadrupole mass filter and ion counting 
was computer controlled using the commercial Extrel Merlin data acquisition and control 
electronics. The system was installed in a vacuum chamber and pumped by a turbomolecular 
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pump, which provided residual gas pressures in the range of 10 7 mbar, even when the 
graphite catcher was heated. 



E. Laser System and Enhancement Cavity 

The excitation and ionization scheme discussed in Section [III CI requires laser systems that 
provide the high power required to saturate the two-photon transition and to efficiently drive 
the ionization steps and at the same time a sufficiently small bandwidth, high stability and 
precise frequency control to reach the required accuracy for the isotope shift measurements. 

For the two-photon transition it is important to have two exactly counterpropagating laser 
beams in the laser ionization region. It is also favourable to have these beams well balanced 
in power to avoid an asymmetry in the background signal due to Doppler-broadened two- 
photon excitation as discussed below. The laser system that fulfilled all these requirements 
was composed of two argon ion laser pumped ring lasers, a titanium:sapphire (Ti:sapphire) 
laser and a dye laser, combined with a Fabry-Perot cavity around the laser interaction 
region to enforce the counterpropagating beams and the high intensities that are required 
for the two-photon transition. Because both laser beams must interact with the atoms 
simultaneously, this solution required a cavity that is kept in resonance with the two laser 
beams having strongly different wavelenghts. High-accuracy frequency determinations were 
achieved by referencing the Ti:sapphire laser to an iodine-stabilized diode laser. 



1. Iodine- Locked Reference Laser 

A stable reference frequency for the isotope shift measurement was realized by an ampli- 
fied diode laser system locked to a hyperfine component in the molecular spectra of 127 l2. To 
obtain a detectable beat frequency between the Ti: sapphire laser and the reference diode laser 
while the Ti:sapphire laser frequency was tuned across the 2s — > 3s transition in lithium, 
the iodine transition had to be within approximately 50 GHz of the two-photon resonance 
frequency. Hence, the X 1 !!* — >■ BO^ i?(114) 11-2 transition in iodine was chosen. The ai 
hyperfine component has a predicted resonance frequency of 407815 138 (30) MHz accord- 



ing to the 'iodine spec 4' program 69L Recent measurements accurately determined the 



frequency as 407815 137.15 (30) MHz 70J. This is very close to the lithium two-photon tran- 
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FIG. 8: (Color online) Diode laser system locked to an iodine line for frequency reference. ASE: 
Amplified spontaneous emission, AOM: Acousto-optical modulator, EOM: Electro-optical modu- 
lator, F: Fiber, FR: Faraday rotator, M: Mirror, P: Polarizer, PD: Photodiode, PID: Proportional- 
integral-differential regulator, PZT: Piezoelectric transducer for tuning the grating of the master 
laser, RAM: Residual amplitude modulation. 

sitions with the isotopes 6,7 ' 8 Li having resonance frequencies below and 9 ' n Li above the iodine 
transition frequency. The largest separation is for 6 Li and is approximately 11.5 GHz. Stabi- 
lization of the diode laser to the iodine transition was achieved using frequency modulation 
(FM) saturation spectroscopy 7l|, with the experimental arrangement shown in Fig. El The 
light produced by the master laser (Toptica, Model PDL 100) was amplified in a broad-area 
diode laser (BAL 740-100-1, Sacher) and separated from the weak amplified spontaneous 
emission (ASE) of the amplifier with a diffraction grating. The spatially-extended amplified 
beam was split at the edge of a mirror into a pump and a probe beam with intensity ratio 
of about 2:1. The pump beam was frequency-shifted (8 MHz) and amplitude- modulated 
(AM) at 6.6 kHz with an acousto-optical modulator (AOM) and then sent through a single- 
mode fiber (Newport S-FS-C) to obtain a good TEMoo mode structure. The probe beam 
was first sent through a mode-cleaning fiber and then frequency-modulated (22 MHz) with 
an electro-optical modulator (EOM) for side-band generation. The signal in FM satura- 
tion spectroscopy after simultaneous interaction of the counterpropagating probe and pump 
beam with the iodine vapour is carried by the amount of amplitude modulation of the probe 
beam at the EOM frequency. Hence, residual amplitude modulation (RAM) of the phase- 
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modulated probe beam introduced by a non-ideal matching of the laser beam polarization 
to the EOM would have caused an offset signal in the detection that shifts the locking point 
away from the resonance center. Therefore the RAM was actively suppressed by separating 
a small part of the laser beam after the EOM and detecting intensity fluctuations at the 

ligh- volt age DC 



72|. 



EOM frequency on a fast photodiode. A feedback loop that regulated a 
offset on the EOM was used to remove the spurious RAM as described in 

The probe and pump beams were superimposed in a counterpropagating geometry with 
perpendicular polarization in an iodine vapor cell (pump beam: 5 mW, probe beam: 3 mW, 
beam diameters 1 mm). The beams were combined and separated with Rochon polarizers. 
After passage through the iodine cell and separation from the pump beam, the probe beam 
was directed onto a photodiode. The photodiode signal was amplified and demodulated with 
a mixer at the 22 MHz EOM frequency. The mixer intermediate frequency (IF) output was 
then fed into a lock-in amplifier for phase-sensitive detection to extract the signal at the 6.6 
kHz frequency that was applied to the AOM driver to amplitude-modulate the pump beam. 

A typical spectrum of the iodine transition obtained by scanning the unlocked diode laser 
across the hyperfine resonances is shown in Fig. [9j To lock the laser to one of the isolated 
hyperfine components, it was coarsly tuned to the frequency of the I2 &i transition guided by 
a wavemeter. The output of the lock-in amplifier was used as input signal for a PID regulator 
to produce a servo signal for the diode laser stabilization. The ai hyperfine component of the 
transition was chosen as the reference transition since it is clearly separated from the other 
hyperfine components and could be easily identified to facilitate reliable relocking during 
the measurements. 

To populate sufficiently the v — 11 vibrational level of the I2 electronic ground state, the 
iodine cell is heated in an oven to 600°C. The iodine reservoir is kept outside the oven in 
a cold finger at a fixed temperature of 29°C to control the vapor pressure inside the cell. 
Pressure broadening of the lines was only observed at cold-finger temperatures well above 
30°C, the pressure-dependent shift of a few kHz/Pa 70[ is not relevant for the accuracies 
targeted in the measurements reported here. 
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FIG. 9: (Color online) Hyperfine structure spectrum of the X 1 ^^" — > BO^ i?(114) 11-2 reference 
iodine line measured with the set-up shown in Fig. [8j 

2. The Titanium: Sapphire Laser 

The 735 nm light for the lithium 2s — > 3s two-photon transition was produced by a 
Coherent 899-21 titanium:sapphire (Thsapphire) ring laser pumped with 15 W from an 
argon ion laser (multi-line visible). The Thsapphire laser provided up to 1 W of single- 
frequency output with a typical linewidth of approximately 1 MHz. Short-term frequency 
fluctuations were suppressed by locking the laser to the external Fabry-Perot cavity that is 
part of the Coherent 899-21. Medium (> 30 ms) and long-term stabilization was achieved 
with a frequency-offset lock of the Ti: sapphire laser relative to the iodine stabilized diode 
laser. About 1 mW of light from the reference diode laser was superimposed with a few mW 
from the Thsapphire laser by coupling each into the two arms of a fiber-optical beamsplitter. 
One of the output arms with the mixed laser beams was connected to the fiber-optical input 
of a fast photodiode (New Focus, Model 1434, 25 GHz), while the output of the second arm 
was collimated and mode matched into a 300 MHz Fabry-Perot interferometer for spectral 
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FIG. 10: (Color online) Optical setup of the Pound-Drever-Hall stabilisation of the dye laser 
(coming from top) and the Ti:sapphire laser (coming from the bottom) and injection of the laser 
beams into the experimental setup for resonance ionization spectroscopy of lithium isotopes. DBS: 
Dichroic beamsplitter mirrors, EOM: Electro-optical modulators for side-band generation, FR: 
Faraday rotators for separation of the reflected laser beams and detection with fast photodiodes, 
L: Mode-matching lenses, A/4 and A/2: Waveplates for polarisation adaption, P: Polarizers, PZT: 
Piezoelectric transducer. 



analysis of the diode and the Thsapphire laser beams. The radio frequency (RF) output 
signal of the fast photodiode, i.e., the beat signal between the Thsapphire laser and the diode 
laser, was amplified with an ultra-wideband amplifier (MITEQ Model JS4-001020000-30-5A) 
and then divided with a two-way splitter. One part was guided to a microwave frequency 
counter (Hewlett Packard, Model 5350B), while the second part was used for frequency offset 
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locking of the Tksapphire laser. To obtain the respective servo signal, the beat frequency 
was mixed with the RF output of a synthesized sweeper (Hewlett Packard, Model 83752A) 
and a frequency difference of exactly 160 MHz between the sweeper and the beat frequency 
was maintained. This was accomplished with a frequency discriminator (MiTeq, Model 
FD-2PZ-160/10PC) that provided an output voltage proportional to the deviation of the 
input frequency from 160 MHz. Its output was fed into a PID regulator and the correction 
signal applied to the external scan input of the Coherent Tksapphire laser. The bandwith of 
the feedback loop was approximately 30 Hz. By changing the frequency of the synthesized 
sweeper, the Ti:sapphire laser could be set and stabilized to any arbitrary frequency around 
the iodine reference line within the bandwidth of the fast photodiode (±25 GHz). An upper 
limit for the laser linewidth could be obtained from the frequency spectrum of the beat 
signal. At GSI it was usually on the order of 1 MHz jl|, while at TRIUMF a slightly larger 
linewidth of 2 to 3 MHz was observed. This increase was mainly due to the frequency jitter 
of the diode laser, caused by the acoustic noise of a CAMAC crate located nearby. However, 
with about 1 s integration time, the average frequency measured by the frequency counter 
was stable to within a few 10 kHz. 

The Ti:sapphire laser light was transported from the laser laboratory to the experimental 
hall with a 25 meter long photonic-crystal fiber (Crystal Fibers, LMA-020) that can trans- 
mit high cw powers without significant losses and without indication of nonlinear Brioullin 
scattering due to its large mode area jisl . Typical transmission through the fiber was about 
80% for 1 W input Ti:sapphire laser power. 

3. Light Enhancement Cavity 

High intensities are required to approach saturation for the 2s — > 3s two-photon transition 
and hence efficient detection. One approach that is often used in Doppler-free two-photon 
spectroscopy is strong focussing of two counterpropagating laser beams to reach saturation 
intensity. However, this has the disadvantage of poor spatial overlap between the two laser 
beams and the atomic beam released from the graphite catcher. Hence, we developed an 
optical enhancement cavity to both ensure collinearity of the counterpropagating beams and 
to obtain higher intensities. A critical point is that both laser beams, one for the two-photon 
resonant excitation and the other for the 2p — > 3d excitation, i.e., the Ti:sapphire and the 
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dye laser beams, must be coupled to the cavity simultaneously. In designing this cavity one 
has also to consider that the spatial profile of the resonator mode should provide a relatively 
large focus to ensure sufficient spatial overlap with the atomic beam. 

Figure [TU] shows the cavity together with the optical setup for stabilization. The symmet- 
rical cavity has a length of approximately 30 cm and mirror curvature radius of 50 cm. The 
diameter of the TEMoo mode in the focus is therefore relatively large and is approximately 
500 /im. The cavity is placed completely inside the vacuum chamber with two mirror holders 
mounted on a vertically oriented baseplate. The input coupling mirror above the ionization 
region of the QMF has a transmittance of 98% while the high reflector mounted below has 
reflectivities of R > 99.93% for both wavelengths (735 nm, 610 nm). The high reflector is 
fixed to a piezoelectric transducer (PZT) for fine tuning of the cavity length. 

The cavity length is actively stabilized to the Ti:sapphire laser frequency using Pound- 
Drever-Hall locking 0, 0]. Collimating optics for the fiber, lenses for spatial mode- 
matching and sideband generation with an EOM required for the locking scheme is mounted 
on a breadboard on top of the vacuum chamber. The Ti:sapphire light is first focused 
through an EOM operated at 12 MHz. Similar to the diode laser modulation discussed 
in Section IIII E 11 an active feedback loop suppresses residual amplitude modulation. Af- 
ter passing through an initial polarizer and Faraday rotator, the main Ti:sapphire beam 
is superimposed on the dye laser beam (see below) using a dichroic beamsplitter that is 
highly reflective for 735 nm at 45° incidence and antireflection coated for 610 nm. Behind 
the dichroic beamsplitter, a broadband mirror directs the combined laser beams vertically 
through an antireflection-coated viewport into the vacuum chamber and to the enhancement 
cavity. Light reflected by the input coupler of the cavity is separated at the input polarizer 
after returning through the Faraday rotator and detected with a fast Si-PIN photodiode. 
Depending on the resonance condition of the cavity, the two sidebands created by the EOM 
modulation exhibit different phase shifts and a dispersionlike signal is obtained after de- 
modulation of the photodiode signal at the EOM frequency. This signal is used to generate 
a PID-regulated servo signal (0-500 V), which is applied to the PZT to change the cavity 
length. This stabilizes the cavity length to the Ti: sapphire frequency and tracks it while 
tuning across the resonance transitions of the different lithium isotopes. The effective tuning 
range of the cavity for PZT voltage variations of 500 V is about 600 MHz and thus suffi- 
cient to cover the hyperfine structure of each lithium isotope (typically 300 MHz) without 
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relocking the servo loop. 

Such an optical resonator is extremely sensitive to vibration because the distance of the 
mirror positions must be stabilized within a small fraction of the laser wavelength, i. e. less 
than 10 nm. To decouple the vacuum chamber with the optical cavity from vibrations 
caused by vacuum pumps and other mechanical devices, the two cross pieces housing the 
experiment (Fig. [T0|) are connected to the beam line only through a flexible metal bellow, 
while the roughing pump is isolated from the chamber turbo pump by flexible tube couplings 
to and from a massive ballast chamber mounted on the floor. 

Fluctuating light intensity inside the cavity has a strong influence on the observed line- 
shape for two reasons: The quadratic dependence of the excitation efficiency increases the 
sensitivity of the signal intensity to vibrationally induced power changes. Even more im- 
portant, changing intensities cause fluctuating ac Stark shifts of the transition frequencies. 
These are discussed in more detail below. Therefore it is important to monitor the light 
power contained within the resonant cavity. To do so, the light transmitted through the 
high reflector, which is about 0.05% and 0.07% of the power inside the cavity for 735 nm 
and 610 nm, respectively, leaves the vacuum chamber through a second viewport. A sec- 
ond dichroic beamsplitter separates the two wavelengths and send them to two photodiodes 
where their intensity is recorded. From this averaged signal a power enhancement of about 
80 to 100 was determined for the cavity. 

4- Dye Laser 

Laser light at 610 nm is required to drive the 2p 2 Pj — > 3d 2 Dj/ transitions in lithium. 
This was produced by a Coherent 699-21 dye laser operated with a Rhodamine 6G solution 
in ethylene glycol and pumped with 6 W of multiline visible output of a second argon ion 
laser. The laser light was transported from the laser laboratory to the experimental hall 
with a second large-mode-area photonic-crystal fiber (LMA-020). In this case, transmission 
efficiencies up to about 70% for 600 mW input dye laser power were achieved. This light had 
also to be resonantly coupled into the optical cavity together with the Ti: sapphire laser light. 
Hence, the dye laser was locked to a longitudinal cavity mode that was as close as possible 
to the 2p 2 P — > 3d 2 D transition frequency. This was achieved using a second Pound-Drever- 
Hall servo loop depicted in Fig. [TUJ but the servo signal was applied to the external-scan 
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input of the Coherent 699 dye laser controller (rather than to the PZT controlling the length 
of the enhancement cavity). This approach does not allow tuning the dye laser exactly to the 
2p — > 3d resonance transitions, but the high laser intensities caused strong power broadening 
of the allowed dipole transition, as will be discussed in the next Section. With typical power 
broadened linewidths of 7 GHz (FWHM) and a free spectral range of the optical cavity of 
500 MHz it was always possible to find a locking point that provided maximum excitation 
efficiency and kept the dye laser frequency in a region of constant excitation efficiency. 

F. Data Acquisition 

Data acquisition (DAQ) of the experiment was based on the Multi-Branch System (MBS) 
developed at GSI. A CAMAC-GPIB controller allowed communication with the RF synthe- 
sizer and the microwave counter. Digital-to-analog, analog-to-digital converters and scalers 
were directly implimented in the CAMAC crate. The MBS system worked stand-alone. It 
read all data from a CAMAC crate and wrote it on a local tape drive. No other computers 
were necessary to take data and to store them. For analysis, the data of each scan was 
transferred via a TCP/IP connection into the data analysis package Origin 7.0, running on 
a PC. Routines written with C++ in the Origin environment were finally used to analyze the 
data. Nonlinear least-square fits of lineshapes were performed with a Levenberg-Marquardt 
algorithm adapted from Numerical Recipes [3]. 

IV. RESULTS AND DISCUSSION 

We present and sumarize results and data from the preparatory experimental phase at 
GSI and TRIUMF and the four beamtimes that were performed at these facilities: an on-line 
beamtime at GSI in December 2003, an off-line beamtime at TRIUMF in June 2004, and 
two on-line beamtimes at TRIUMF in September and October 2004. 

A. Lineshapes 

Determination of the isotope shift on the level of 100 kHz or better requires a detailed 
understanding and description of the resonance lineshape and all factors affecting it. Hence, 
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this Section starts with an analysis of the observed lineshapes and the influence of the ac 
Stark effect in the 2s 2 Si/2 —> 3s 2 Si/2 as well as in the 2p 2 Pj — > 3d 2 Dj/ transitions. 
Afterwards, the measurements of the hyperfine structure (HFS) and isotope shift (IS) of the 
stable and short-lived isotopes are discussed. 

1. 2s — > 3s Two-Photon Transition 

Resonance profiles of the 2s — > 3s transition were recorded in the following way: First, 
the Ti:sapphire laser was tuned to a frequency below the 2s — > 3s resonance of the respective 
isotope and frequency-offset locked to the iodine-stabilized diode laser. Then, the enhance- 
ment cavity was locked to the Ti:sapphire laser and finally the dye laser was locked to the 
longitudinal mode of the cavity closest to the 2p — > 3d transition frequency of the respec- 
tive isotope. The Ti:sapphire laser was scanned by slowly varying the RF frequency for the 
frequency-offset lock in steps of 1 MHz. Maximum scan ranges were about 600 MHz, limited 
by the voltage range of the high-voltage power supply for the piezoelectric transducer (PZT) 
at the enhancement cavity. Since the HFS of all lithium isotopes in the two-photon tran- 
sition is smaller than this scanning range the complete resonance structure can be covered 
with a single scan. 

Figure [TT] shows an overall resonance profile obtained for 7 Li. The observed count rate at 
the channeltron-type detector is plotted as a function of the measured beat frequency be- 
tween the iodine-locked diode laser and the Ti:sapphire laser. Here, a large tuning range was 
used to record also the far-reaching wings of the resonance profile. Therefore the complete 
scanning range of about 2 GHz was divided into scans of approximately 400 MHz range in 
order to stay inside the tuning range of the cavity PZT and to limit the frequency variation 
of the dye laser (see below). After each 400 MHz scan, the cavity and the dye lasers were 
relocked to keep the dye laser frequency as close to the 2p — > 3d resonance frequency as 
possible. 

The profile exhibits two narrow Doppler-free components that are labeled with their F 
quantum numbers. According to the selection rules for a two-photon transition between s 
states, only hyperfine transitions with AF = are allowed. The wide Doppler-broadened 
background with a relative amplitude of about 1% of the peak intensity and a Gaussian width 
of approximately 1.7 GHz is caused by the absorption of two photons from copropagating 
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FIG. 11: (Color online) Wide-range line profile of the 7 Li 2s — > 3s two-photon transition with 
narrow Doppler-free hyperfine transitions and the Doppler-broadened background as observed off- 
line at GSI. On the x-axis, the beat frequency between the iodine-stabilized diode laser and the 
Ti:sapphire laser is given. The solid (red) line is a fit of two Voigt profiles including a wide 
Doppler background pedestal. The dashed (green) line shows the individual contributions of the 
two hyperfine components. The regions marked by vertical lines are the scanning regions in 7 Li 
on-line experiments. The inset shows schematically the hyperfine splitting in the 2s and 3s states 
of the I = 3/2 isotopes ,9,11 Li and the allowed (AF = 0) two-photon transitions. 



beams 



77l | . To account for this process, each of the two peaks is fitted by applying a 



Levenberg-Marquardt minimization procedure with a Voigt plus a background Gaussian 
function represented by the dashed (green) lines. Both HFS components are constrained 
to have identical width parameters for the Voigt as well as the Gaussian lineshape of the 
background. The Voigt profile shows a Lorentzian linewidth I\ of 4.5 MHz which is already 
slightly larger than the natural linewidth of 2.6 MHz, attributed to saturation broadening, 
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TABLE VII: Estimated partial efficiencies and expected overall efficiency for the detection of 
lithium ions compared with that experimentally observed. 


Release cffirioncv of the catcher foil 


50 % 


1 lifpr an nr»'l" , ixrooTi mQpr npaniQ ann afnynip ovarii (Hi mm II ^ lmim 1 




T^'XTit^tion pfrir'iP'nr'v tot flip — V 3<3 two-'nlint on tv^'nmtin'n 


25 % 


Ionization efficiency via 2p — )• 3d — > continuum 


3 % 


Signal reduction by HFS splitting 


62 % 


Transmission of the quadrupole mass filter 


90 % 


Quantum efficiency of the detector 


80 % 


Expected overall efficiency on resonance 


2.5 • I0~ 4 


Experimental overall efficiency at GSI 
Experimental overall efficiency at TRIUMF 


2 • 10~ 4 
5 • 10~ 4 



and a Gaussian linewidth Tq of 1 MHz, fitting well to the observed laser linewidth. Here, 
all width parameters refer to the Tbsapphire laser frequency scale and have to be doubled 
to obtain the value for the 2s — > 3s two-photon transition. The solid (red) line is the overall 
fitting function which shows an excellent agreement with the experimental data points over 
more than 2 GHz frequency range and about three orders of magnitude in signal intensity. 

For precise frequency determination, only the regions around the two Doppler-free peaks 
are important and scans of the radioactive isotopes were thus usually performed by scanning 
about ±30 or ±15 MHz around the resonance centers, as indicated in Fig. [TTJ skipping the 
intermediate part in fast steps without data taking. Such spectra of 6 Li and 7 Li, taken at 
the GSI on-line mass separator and the Off-Line Isotope Separator (OLIS) at TRIUMF, 
respectively, are depicted in Fig. [12] The incoming ion yield was about 3 • 10 5 ions/s for 6 Li 
and 2 • 10 8 ions/s for 7 Li, while approximately 80 and 10 5 ions/s were obtained in resonance 
on the strongest hyperfme transition. This corresponds to overall efficiencies of ~ 2 ■ 10~ 4 at 
GSI and 5-10 -4 at TRIUMF. The observed efficiencies agree quite well with those calculated 
and estimated during the design of the setup as listed in Table I VI II 

Before fitting the recorded spectra, the observed number of ion events (iVR aw ) was nor- 
malized for each channel with the Ti: sapphire laser power (PTi:sapphire) that was recorded 
with the photodiode located behind the enhancement cavity (see Fig. ITUl) while the atoms 
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were irradiated with the constant laser frequency: 



AW = AW ■ ( ^ Ti:sapphire ) . (64) 

V V Ti:sapphirc/ / 

Here (-PTi:sa PP hirc) is the average Ti:sapphire laser power while recording the complete spec- 
trum. The uncertainty of A^orm was calculated from the statistical uncertainty and from 
the laser power fluctuations using Gaussian error propagation. To confirm the normaliza- 
tion function, a x 2 optimization was performed using different normalization exponents and 
checking for the lowest x 2 value in the subsequent lineshape fitting. The optimum was found 
to scatter between 1.7 and 2.8, and hence the theoretically expected factor of 2 seems to be 
well justified. 

The red fit curves are Voigt profiles including a Doppler background with a width fixed 
to 1.7 GHz as obtained from Fig. [TTJ The Lorentzian and Gaussian linewidths of the narrow 
Voigt profile as well as the relative intensity of the Doppler background are constrained to 
be equal for both hyperfine structure components. The width of the Gaussian pedestal was 
varied within a reasonable range (~ 1 - 3 GHz) and does not show considerable influence 
on the fitted peak centers. The fitting curves of the individual peaks are indicated by the 
green dashed lines. This procedure was used for all isotopes. 

In the lower part of Fig. [121 the residuals of the fit are plotted. A slight systematic 
asymmetry is visible in the 7 Li spectrum while it is much less pronounced in the residuals of 
the 6 Li resonance fit. In later experiments performed off-line at GSI, a much stronger and 
clearly visib 
publication 

the laser beam. Atoms passing through the laser beam in the interaction region experience 
an intensity-dependent ac Stark shift as will be discussed in the next section. Therefore, 
the exact atomic resonance frequency depends on the position within the laser beam where 
the excitation and ionization occur. This line profile distortion is similar for all isotopes 
under identical experimental conditions. Therefore, the line shape calculations described in 
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e asymmetry of the peaks was observed which is discussed in detail in a recent 
2^| . There it is shown that the asymmetry is caused by the Gaussian profile of 



27| show that the influence of this asymmetry on the extracted isotope shift is very small. 



However, its contribution is significant if the total transition frequency is to be extracted. 



Then a total correction by about 160 kHz is required [27]. 

To calculate the isotope shift, the resonance positions of the individual hyperfine compo- 
nents obtained from the fit must be converted into center of gravity (e.g.) frequencies. Since 
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FIG. 12: (Color online) 2s — > 3s two-photon resonances for (a) 6 Li obtained at GSI and for (b) Li 
recorded at TRIUMF. The x-axis break indicates the intermediate region that was skipped. The 
solid (red) line is the common fit function while the dashed (green) lines are the contributions of 
the individual hyperfine components. The residuals (differences between the experimental points 
and the common fit function) are plotted at the bottom. 



the nuclear quadrupole moment does not affect the J = 1/2 states in the transition, only 
the magnetic hyperfine interaction must be considered. The energy shift of the hyperfine 
state with angular momentum F relative to the J-level energy is in first order given by 



£hfs =2 Cf = ^ + X) " J( J + X) _ J(J + X) ] 



(65) 



with the Casimir factor Cp and the magnetic dipole hyperfine constant A. In first-order 
perturbation theory, the hyperfine structure eg coincides with the unperturbed J-level energy 
and can be calculated from the two hyperfine resonances in the 2s — > 3s transition of lithium 
according to 



Cp vpi — Cpi fp 



(66) 



Cp — Cpi 

where Vp is the transition frequency of the F — > F transition. For 7 ' 9 ' n Li with nuclear spin 
/ = 3/2, this leads for example to is cg = |z/ 2 + f^i- 

Figure [13] shows the Lorentzian linewidth of the two-photon transition as a function of 
the Thsapphire laser intensity inside the cavity (/Tiisapphirc) • The intensity was calculated 
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from the power transmitted through the high reflector of the enhancement cavity with a 
measured cavity mirror transmission of 0.05(1)% and a calculated cavity mode diameter of 
0.46 mm. The solid red line is a fit of the function for the power-broadened linewidth of a 
two-photon transition 

Tl = T \Jl + ((/Ti:sapphirc)//Sat) 2 (67) 

to the data points, where T = 3.08(2) MHz is the extrapolated natural linewidth for 
vanishing intensity of the Ti:sapphire laser and /sat = 148(2) W/mm 2 is the saturation 
intensity. These results are in good agreement with values obtained in a previous beamtime 
at GSI, where r = 3.2(1) MHz and ig at = 167(6) W/mm 2 were obtained. The uncertainty of 
the absolute value of /sat is solely the fitting uncertainty and includes neither the uncertainty 
in the mirror transmission and the transmission through the entrance window 4 nor the 
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FIG. 13: (Color online) Lorentzian linewidth of 6 Li resonances as a function of the Ti:sapphire 
laser intensity as obtained in off-line experiments at TRIUMF. The solid line is a fit by Tl = 
royl + (7/^Sat) 2 ) describing the saturation for a two-photon transition. The Ti:sapphire laser 
intensity is calculated based on the laser power behind the high reflector of the enhancement 
cavity as described in the text. Please note that the linewidth refers to the Ti:sapphire laser 
frequency and must be doubled for the full 2s —> 3s transition. 

4 The window transmission shows an etalon effect discussed below and changes slightly as a function of the 
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uncertainty of the effective diameter of the laser focus. Hence an additional uncertainty 
of about 20% should be added if these values are to be compared with calculations or 
other measurements. The small deviation of the fit value T from the natural linewidth 
Tonat = 2.65 MHz is either caused by laser intensity fluctuations inside the cavity, resulting 
in varying ac Stark shifts for the atoms of the ensemble, or an artefact from the fit program. 

2. 2p — > 3d Resonance 

The lineshape of the two-photon transition is distorted if the ionization efficiency along 
the 2p — > 3d — > Li + ionization path changes during the scan. This could happen because the 
change in the Thsapphire laser frequency during a scan of the 2s — > 3s transition induces 
also a change of the dye laser frequency. This cannot be avoided since the dye laser is locked 
to the enhancement cavity and the resonator length must be changed when scanning the 
Ti:sapphire laser frequency. Hence, the dye laser frequency is changed typically by a few 
hundred MHz during the scan. Helpful is here the strong power broadening of the 2p — > 3d 
transition. The high intensity of the 610 nm light in the cavity focus considerably broadens 
the linewidth of the 2p — > 3d transition. This broadening was studied to find maximum 
excitation efficiency and to ensure a constant ionization efficiency along the whole scan. 
The results for different laser powers as they were available at TRIUMF and GSI are shown 
in Fig. HH The main reason for the different power levels is the fiber transport from the 
laser laboratory to the experimental hall. While a 50 m long standard single-mode fiber was 
used for the dye laser at GSI, a large mode area (LMA) photonic crystal fiber was applied at 
TRIUMF with a transport distance of 25 m. Hence 20 mW of dye laser light were coupled 
to the enhancement cavity at GSI while 80 mW were available at TRIUMF. 

The resonance profiles in Fig. [T4"lwere obtained in the following way: First the QMF was 
set to detect photo ions of the respective isotope and the Ti:sapphire laser frequency was 
fixed at the resonance of the strongest hyperfine transition and operated at about 25% of the 
maximum achievable power. Then, the frequency of the dye laser was set to a value clearly 
below resonance and changed manually until a longitudinal mode of the enhancement cavity 
was reached. There, the dye laser frequency was locked to the cavity and resonant laser ions 

laser wavelength. 
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FIG. 14: (Color online) Strongly saturated resonance profiles of the 2p — > 3d transition in Li as 
obtained with 20 mW laser power entering the resonator at GSI (a) and with 80 mW at TRIUMF 
(b). Frequencies were measured with a Fizeau-type interferometer and have uncertainties of about 
3 • 10~ 3 cm -1 . Circles indicate the experimentally obtained signal and the solid line is a fit of 
two Lorentzian line profiles with identical widths to the data. This curve is only to guide the eye 
since the strong saturation and the unresolved 2 D 3/ / 2j 5/2 levels lead to a distorted line shape. The 
flat-top regions used for spectroscopy in the 2s — > 3s transition are indicated by vertical lines. The 
inset in (b) is the result of another scan with higher statistics around the flat-top region. 

were detected for a period of 3 s before the dye laser was taken out of lock and its frequency 
changed until the next longitudinal cavity mode was reached. The free spectral range of the 
cavity was approximately 500 MHz which determined the step size. The intensity of the 
incoming Li + ion beam as well as of the Thsapphire laser was sufficiently constant during 
the measurements and a normalization of the count rate therefore not required. 

In Fig.[14]the detected ion count rates (circles) are plotted as functions of the dye laser 
frequency that was determined with a Fizeau-type wavemeter. The solid line is a fit of 
two Lorentzian profiles with equal widths to the data and serves only to guide the eye. 
The profile (a) obtained at GSI with lower dye laser power shows a flat-top region with 
an extension of about 1.5 GHz. That is sufficient to ensure constant ionization efficiency 
since the scan width of the Ti: sapphire laser was less than 500 MHz for all isotopes. In 
(b), recorded with the higher dye laser power that was available at TRIUMF, the flat 
region extends to approximately 0.1 cm -1 (3 GHz) width, ideal for looking the dye laser 
during a scan. The region of the first resonance agrees roughly with the reported transition 
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TABLE VIII: Mass shift (^p-adMS an< ^ resonance frequencies i>2p-Ad of the 2p 2 P&2. ~ * 3d 2 D 5 / 2 
transition in lithium isotopes based on the 6,7 Li transition frequencies reported in 78] 



Isotope 


6 Li 


7 Li 


8 Li 


9 Li 


n Li 


Su 2p-3d,MS ( MHz ) 
V2p-3d 




16 379.0089 


2 794 
16 379.1021 


4 901 
16 379.1724 


6 535 
16 379.2269 


8 918 
16 379.3064 



frequencies of 16 379.0661 cm" 1 and 16 379.1021 cm" 1 [78| for the 2p 2 P 3/2 3d 2 D 3/2 
and 2p 2 P 3 / 2 — > 3d 2 D 5 / 2 resonances in 7 Li, respectively. Obviously, the dye laser power 
at TRIUMF was sufficient to even excite the tails of the 2p 2 Pi/ 2 — > 3d 2 D 3 / 2 transition at 



16 379.4014 cm -1 [78[. Contrary, in the measurements at GSI (Fig. [T47 a)). where less dye 
laser power was available, the 2p 2 Pi/ 2 — > 3d 2 D 3 / 2 transition is still clearly resolved. Thus, 
we can expect the overall efficiency at TRIUMF to be about 30% larger than that obtained 
at GSI since the dye laser will also excite and ionize the fraction of atoms decaying from the 
3s 2 Si/ 2 to the 2p 2 Pi/ 2 state. This is in accordance with the higher efficiencies obtained at 
TRIUMF as reported above. 

The frequency of the 2p 2 P 3 / 2 —> 3d 2 D 5 / 2 resonances for the unstable isotopes were 



calculated from the isotope shift formula 



8v 



AA' 



K 



M A M A , 



MS ' 



(68) 



with K MS = 117.8 
isotope shift data 



M A - My 

2 GHz-amu being the mass shift coefficient obtained from the known 6 ' 7 Li 



and neglecting all field shift contributions. The calculated resonance 



frequencies that were used for setting the dye laser frequency are listed in Table IVIIII 



B. Ac Stark Shift 



Since atoms experience relatively strong fields when crossing the laser beam focus in 
the resonator, the atomic level energies will be altered by ac Stark shifts. To correct for 
these shifts, spectra were recorded at different light powers. To distinguish between the 
effects caused by the Ti: sapphire and the dye laser light, the power of one of the lasers 
was varied while keeping the other laser intensity constant during the measurements. The 
observed eg frequencies were plotted against the power level of the respective laser detected 
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FIG. 15: (Color online) Ac Stark shift in the 2s — > 3s transition of 6 Li caused by (a) the Ti:sapphire 
laser light at 735 nm and (b) the dye laser light at 610 nm. In (a) the power of the dye laser was 
fixed to 50 mW in front of the cavity and in (b) the power of the Ti:sapphire laser was fixed to 
330 mW during the measurements. The solid lines are linear regresion lines while the dashed lines 
indicate the 95% confidence bands of the fit. 

on the photodiodes behind the cavity high-reflector. These values were not converted into 
intensities at the laser focus since this conversion includes a large uncertainty due to the 
mirror and window transmission functions and the effective laser diameter inside the focus. 
Figure IToT a) shows the 6 Li eg frequency as a function of the Thsapphire (a) and the dye 
laser power (b). In both cases a linear relation is observed as is expected for an off-resonance 
ac Stark shift. Hence, the linear function 

Vc^bo + hiP) (69) 
is fitted to the data points from which the slopes 

6 Ti: S a PP hirc = _ a0549 ( 14 ) MHz/mW and (70) 

bi yc = 0.115(14) MHz/mW (71) 

are obtained. In both fits of Fig. [T^l the 95% confidence bands are indicated. The dye and 
the Ti:sapphire laser intensities influence the level energies in opposite directions, resulting 
in the different signs of the slope. Even though the sensitivity on the dye laser intensity is 
stronger, the absolute shift that it causes is much smaller due to the relatively low power 
levels compared with the Thsapphire laser. Varying from highest to lowest applied power, 
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the Ti: sapphire laser light shifts the 2s — > 3s frequency by about 1.5 MHz, the dye laser light 
only by 150 kHz. Hence, it was concluded that it is sufficient to keep the dye laser power 
as constant as possible during all measurements whereas such power series measurements 
were performed with the Tksapphire laser for all isotopes excluding n Li. For the latter 
isotope the production rates were not high enough to produce sufficient statistics at low 
laser power. A remaining small influence of the ac Stark shift, caused by power variations of 
the dye laser while taking the Ti:sapphire laser power-series, was included in the uncertainty 
of the isotope shifts as discussed below. 



C. Hyperfine Structure and Isotope Shifts 

A large number of measurements on the stable isotopes was performed during the prepa- 
ration of the experiments and the four beamtimes that were performed at GSI and TRIUMF. 
Spectra like those shown in Fig.lT2lwere analyzed by fitting the lineshapes as discussed above 
and determining the e.g. frequency z/ cg and the hyperfine splitting Az/^ FS from the peak po- 
sitions. Spectra of the radioactive isotopes will be presented in the following Sections. 



1. Hyperfine Structure 

The hyperfine structure splitting in an atomic transition is determined in first-order 
perturbation theory by the following nuclear ground state properties: spin J, magnetic 
moment /i, and the spectroscopic electric quadrupole moment Q s . Values of those are listed 
for all lithium isotopes in Table IIXI Additionally, the magnetic hyperfine constant A 2s in 
the 2s 2 Si/2 electronic ground state, which is directly proportional to the magnetic moment, 
is included in the table. For the stable isotopes they are well known from atomic beam 
magnetic resonance measurements [3| and for the short-lived isotopes they are calculated 
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22]. 



based on the ratio of the magnetic moments taken from 

The spectra of the two-photon transitions recorded here, provides only information about 
the hyperfine splitting Ai/ HFS between the F — 1+1/2 -+ F — 1+1/2 and the F — 1—1/2 ->■ 
F = I — 1/2 transitions. Separate determinations of the hyperfine splitting in the ground 
and the excited state are therefore not possible. However, using the information listed in 
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TABLE IX: Nuclear spins and parity I 77 , experimental dipole nr and spectroscopic quadrupole 
moments Q s of the investigated lithium isotopes from 
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221 ] . A 2s values are taken from 



79] 



for the stable isotopes and are calculated based on the ratio of the magnetic moments for the 
short-lived ones. 



Property 



*Li 



7 Li 



*Li 



9 Li 



ii 



Li 



r 

V-i (mat) 
Qs (mb) 
A 2s (MHz) 



1+ 

0.8220473(6) 

-0.806(6) 

152.136840(3) 



3/2" 

3.2564268(17) 

-40.0(3) 

401.752044(3) 



2+ 

1.653560(18) 

+31.4(2) 

153.017(3) 



3/2" 

3.43678(6) 

-30.6(2) 

424.031(7) 



3/2" 
3.6712(3) 
(-)33.3(5) 
452.954(37) 



Table HX| the hyperfine splitting in the excited 3s 2 Si/ 2 level can be extracted according to 

A..HFS 

A 3s ( A U) = A 2s ( A U) - (72) 

and the results will be used to check for an indication of hyperfine anomaly which might be 
expected for such an extended nucleus like n Li. 

Results for A^ HFS are plotted in Fig. [TBI Here, each data point is the average of typically 
10 to 20 measurements with the standard error of the mean (SEM) as statistical uncertainty. 
The rightmost data point is the overall weighted average, where the standard deviation of all 
results is represented with the dashed lines and used as the final uncertainty. Additionally, 
X 2 -fitting of a constant to the data points was performed. It resulted in reduced-^ 2 values 
between 0.5 and 2.0 for all isotopes, indicating that the SEM of such a set of measurements 
is a realistic estimation of the corresponding statistical uncertainty. Hence, for n Li, for 
which in total 24 spectra were obtained (see below), the corresponding SEM is used as 
uncertainty of the weighted average of all individual measurements. The average values and 
uncertainties obtained in this way are listed in the first line of Table Kl 

The A factors of the 3s state were calculated using Eq. ( 1721) and are listed in in the last 
line of Table |2D The current result for ^4 3s ( 7 Li) disagrees with the previous experimental 
determination by Stark spectroscopy 80J, but is in excellent agreement with a previous laser 



spectroscopic measurement on an atomic beam 25j and with theoretical values, obtained by 
both Hylleraas variational 8l| and multiconfiguration Hartree-Fock methods £ 



. After the 

measurements at TRIUMF, frequency-comb based measurements of the 2s — > 3s transition 
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FIG. 16: (Color online) Results of all hyperfine splitting measurements in the 2s 2 Si/2 — > 3s 2 Si/2 
transition of Li isotopes as obtained during the off-line and on-line beamtimes at GSI and TRIUMF. 
For 6 ' 7 Li, additional off-line measurements were performed at GSI after the TRIUMF measurements 
using a frequency comb for laser frequency stabilization j^]. These results are shown as the 
rightmost GSI datapoint and are in good agreement with previous data. 

Plotted are the measured two-photon frequency intervals Au^ FS (twice the Ti:sapphire laser 
frequency difference) between F — >■ F transitions for F = I — 1/2 and F = I + 1/2. Each data 
point (circle) is the average of 10-20 measurements. The rightmost data point (full circle) 
represents the weighted average. For their uncertainties see caption of Table [ID 



53 



TABLE X: 2s — > 3s hyperfine splitting Ai/ HFS (twice the difference in Ti:sapphire laser resonance 
frequencies between the F = I + 1/2 -)• F = I + 1/2 and the F = I - 1/2 ->■ F = I - 1/2 
transitions) of the lithium isotopes and calculated A 3s constants according to Eq. (|72p . The A 2s 
constants listed in Table IIXI are used. Uncertainties of the average are the geometric sum of 
the statistical uncertainty and a 10 kHz systematic contribution caused by a possible unresolved 
Zeeman effect in the stray magnetic field at the beamline as discussed in Section TlVDI Theoretical 
values Am^Io f° r the hyperfine splitting are obtained from a calculation of the Fermi contact term 
811 ] and experimental nuclear magnetic moments. They include a finite nuclear size (Zemach) 



correction with exponential (e Ar ) and Gaussian (Gauss) electromagnetic distributions 
discussed in the text. 



8l|] as 





6 Li 


7 Li 


8 Li 


9 Li 


n Li 


Al^ p S (Average) 


175.304(21) 


617.301(22) 


293.853(34) 


651.468(29) 


696.09(10) 


A^fs (pnnl) m 
Al/ Theo (Rc oc Gauss) 


175.311(24) 
175.295(17) 
175.305(17) 


617.291(22) 
617.264(60) 
617.289(60) 


293.851(29) 
293.867(23) 


651.470(65) 
651.502(64) 


695.88(9) 
695.92(9) 


A 3s 


35.267(14) 


93.103(11) 


35.476(14) 


98.297(16) 


104.91(6) 



were performed at GSI using the same apparatus but a considerably changed laser system, 
where the Thsapphire laser was locked to a frequency comb 27j]. The factors obtained 
from these measurements are in excellent agreement with the prior ones and are included 
in Fig. [TBJ Also, the ratio A 3s ( 7 Li)/A 3s ( 6 Li) = 2.6407(10) is in good agreement with the 2s 
ground state ratio of 2.6407282(1) from Ref. 83j, indicating no hyperfine anomaly at the 
level of precision available in these experiments. The A 3s factor of the short-lived isotopes 
8 ' 9 Li as obtained from our measurements are also listed. 

According to 8l|], the magnetic dipole coupling constant in a 2 Si/2 state can be written 



as 



A V2 (exp) = 95.410 67(7) 



QeHlft 



exp 



31 



(73) 



in MHz. The experimental Fermi contact term / exp can now be extracted from the known 
magnetic hyperfine constant A 2s . According to theory / exp is connected to the uncorrected 
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Fermi contact term f c by the factorization: 

/exp = — " -CrelCjuCiiCQED/c- (74) 

9e 

Here, f c is the fermi contact term as obtained from nonrelativistic Hylleraas variational 
calculations for a point dipole, a e denotes the anomaly of the electron magnetic moment, 
C re i is a relativistic correction factor, Cm the finite mass correction, Cr the finite size 
correction and Cqed are QED corrections other than the a e . The finite size effect originates 
from the nuclear moment distribution in the nucleus and cannot be evaluated independently 
of a nuclear model. For a one-electron ion, the correction is given by the Zemach correction 

C R = 1- 2Z(R cm )/a (75) 

with the Bohr radius do and the average electromagnetic charge radius (R em ) for the nu- 
cleus obtained by folding the charge and magnetization distributions with sizes R c and i? m , 
respectively. The result for (R em ) depends somewhat on the model chosen for the distribu- 
tions. To calculate the size of possible effects, two models are applied that are discussed also 
in 8l|. In both cases it is assumed that R c = R m and either an exponential distribution 
e -Ar or a Q auss i an distribution is used. This results in different prefactors and therefore 
different values for (R em ) 

35 4 

(-Rem) cxp = ~~~ 7^^ c an d (-Rem) Gauss = -f=R c . (76) 

16v3 V37T 
The charge radii obtained from our measurement and discussed in the next Chapter as 
well as the magnetic moments from literature (Table HXj) allows us to calculate the Zemach 



contributions [81] for all isotopes and the hyperfine splitting Au^^ Q (2s — > 3s). The results 
are included in Table |2D The influence of the two models is small and no model can be 
favored within the accuracy of our results. Even in the case of the halo nucleus n Li where the 
assumption R c = R m is for sure not fulfilled, the values for the theoretical and experimental 
HFS splittings lie almost all in the range of the la uncertainties. 

2. Isotope Shift of ^ 7 Li 

The isotope shift 5v^ between the two stable isotopes was extracted from measurements 
of the dependence of 5vi$ on the Ti: sapphire laser power and by calculating the ac Stark 
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shift regression lines for both isotopes as discussed above. A typical result is shown in 
Fig. [T71 Note that the y-axis frequency scale in the ac Stark shift plots is always given 
in terms of the Ti:sapphire laser frequency. Hence, the isotope shift in the 2s — > 3s two- 
photon transition is twice the distance between the y-axis intercepts of the two ac Stark 
shift regression lines, as indicated in the figure. At this point, the laser power is zero and 
the atomic level energies are unperturbed by the laser light. For both regression lines the 
95% (2a) confidence interval (CFI) was also calculated and is plotted in the graph (dashed 
lines). The statistical uncertainty (la) of the extracted isotope shift was determined from 
the fitting result (lcr confidence limit) and an additional uncertainty of 10 kHz was added to 
account for dye laser power fluctuations during the measurements of the Ti: sapphire-induced 
ac Stark shift. This value was obtained using the b^ ye coefficient according to Eq. ( ITTj) and 
assuming dye laser power fluctuations up to 1/3 of the total power as a conservative estimate. 
Such fluctuations might tilt the ac Stark shift line systematically in a single set. However, 
they will fluctuate statistically between different sets and the further statistical treatment 
is expected to be justified. 

Results of all isotope shift measurements performed at GSI and at TRIUMF are summa- 
rized in Table |XT] and plotted in Fig. [THJ Each entry corresponds to an extrapolation from 
an ac Stark shift measurement as it is shown in Fig. [T71 The isotope shifts obtained in the 
individual sets at TRIUMF were averaged, weighted with the uncertainty, and the weighted 
error is taken as the final statistical uncertainty. Fitting a constant to the data points, as 
shown in Fig. [18] (bottom graph for 6,7 Li), results in a reduced-^ 2 between 0.5 and 1 for all 
isotopes. This x^ed was n °t used to scale the final uncertainty. Excellent agreement between 
the TRIUMF average and the previous measurements at GSI is found for the 6 ' 7 Li isotope 
shift. The combined values from all ToPLiS experiments, including those at GSI, are listed 
in the table and marked 'Average" . 

Since the dye laser power fluctuates statistically, its influence on the line position was 
added to the statistical uncertainty. Other effects may contribute with systematic shifts 
in a preferred direction. The estimated size of these systematic uncertainties is given in 
the caption of Table IXIl It includes contributions from the small lineshape asymmetry 
mentioned before (10 kHz) and possible contributions from unresolved Zeeman shifts caused 
by stray magnetic fields (10 Hz) which will be discussed in Section TlV Dl 

The isotope-dependent velocities cause a differential second-order Doppler shift (SODS) 
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FIG. 17: (Color online) Determination of the isotope shift from the ac Stark shift measure- 
ment. The frequency offset | ^Ti:sapphire — ^dl| between the Ti:sapphire and the iodine- locked diode 
laser at the center-of-gravities of the 6 Li (upper line) and Li (lower line) hyperfine structures is 
plotted as a function of the Ti:sapphire laser power. The x-axis is the average Ti:sapphire laser 
power observed on a photodiode behind the enhancement cavity while recording the resonance. 
This power is proportional to the light intensity at the waist of the resonator mode where the 
atoms cross the laser beams. It is varied between the maximum obtainable value and 25% of the 
maximum value. The dye laser power was kept constant at the maximum obtainable power during 
all measurements. A linear shift of the resonance positions is observed which is caused by the ac 
Stark shift effect. The solid lines are linear fits to the data points and the dashed lines indicate 
the 95% confidence bands. The isotope shift 5zA 7 in this two-photon transition was calculated by 
taking twice the distance between the intercepts with the y-axis as indicated by the arrow. 



57 



between the isotopes 6 Li and A Li with mean square velocities (t>f) and (v\) according to 

5u e,A _ vq {va)-(vI) m) 

°^sods — o „2 y' ') 



Here, we used the value (v\) = 3/2 • 2ksT /Ma for the most probable velocity of the atoms 
in the "quasi" -collimated atomic beam that enters the laser focus. A catcher temperature of 
T « 2000 K is estimated and z/q = 815.61 THz is the 2s — > 3s transition frequency (2 • z/TiSa) 
27(| . Due to the uncertainty of the absolute temperature and the insufficient knowledge of 
the exact velocity distribution, an uncertainty of about 30% is assumed for this correction. 
Absolute shifts are on the order of 30-60 kHz and the differential values are included in 
Table IXII The largest effect appears naturally for the pair 6 Li - 11 Li, which amounts to 
^sods = 34(11) kHz. 

The values reported here were carefully reanalyzed with respect to statistical and system- 
atic uncertainties including effects that were discovered and evaluated I27J after the original 

n 

publication |3[. Hence, the values and uncertainties do slightly deviate from those given 
in our earlier work. However, these changes are small and do not significantly alter the 
determination of nuclear charge radii and their final uncertainty. Note also that all isotope 
shifts listed in Table IXII are given relative to 6 Li, whereas they were earlier given relative 
to 7 Li. The reason for this is twofold: First, the data in Table IXII are the direct results 
since we have used 6 Li as the reference isotope during the measurements. Second, as will be 
described in a following publication, a reevaluated nuclear charge radius of 6 Li appears to 
be more reliable than the previously used reference radius of 7 Li. 

The final value of the ToPLiS measurements for the 6,7 Li isotope shift agrees with the 



result obtained at GSI 



2fl but differs from a previously reported one 25| by about four 
times the combined uncertainties as listed in the last line of Table IXII This is attributed to 
unaccounted systematic errors in the prior interferometric measurements. Possible sources 
of systematic errors are nonlinearities in the scan of the confocal Fabry-Perot interferometer 
(FPI) or changes in the effective length of the interferometer with pointing and wavelength. 
The latter might be caused, for example, by a wavelength dependent phase change of the 
multilayer dielectric mirrors of the FPI. Moreover, the FPI was calibrated at the reference 
lines of Rb at 778 nm and Cs at 851 nm, while the operational wavelength for the measure- 
ment was relatively far away at 735 nm for lithium. A relative error of 2 • 10~ 5 in the length 
of the FPI is sufficient to explain the difference in uncertainty. The data reported here are 
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based on beat frequency determinations and do not rely on any interferometer calibration. 



Please note that the interferometer scale error in |25| is only relevant for the reported isotope 
shift, whereas its contribution to the A 3s factor is only a small fraction of the previously 
stated statistical uncertainty and thus requires no correction. 



3. Extraction of 5 (r 2 ) for 6,7 Li 



We can now use the measured isotope shift between the stable isotopes to extract the 
change in the mean square nuclear charge radius. This is an important test of our tech- 
nique, because the result must be consistent with corresponding values extracted from other 
transitions or by other means, as long as those are model-independent. Therefore, we use 



the calculated mass shift 8v\ 



•J 

MS.Theory 



from Table II III to extract by use of Eq. ( 162]) the change 



in the ms nuclear charge radii and obtain 



K<- 2 > 



2\6,7 



-0.731(22) fm 2 . 



Nuclear charge radii of the stable lithium isotopes were measured bye 



3-861 



86 1 and results are summarized in 
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(78) 

astic electron scat- 
We used the 



tering in the 1960's and 70's 
charge radii _R c ( 7 Li) = 2.39(3) fm 84J and _R c ( 6 Li) = 2.55(5) fm, the average of the values 
reported in 85] and 86|. In Table IXiTl the resulting value for § (r 2 ) 6 ' 7 = —0.79(29) fm 2 is com- 
pared with optical measurements of three transitions in lithium that have been repeatedly 
investigated during the last decades. Besides the 2s — > 3s two-photon transition investigated 
here, isotope shifts were determined in the resonance lines 2s 2 Si/2 — > 2p 2 Pi/2 (Di line) and 
2s 2 S!/ 2 — > 2p 2 P3/2 (D 2 line) in neutral lithium as well as in the ls2s 3 Si — > ls2p 3 Po,i,2 
transitions in Li + . We use the recent mass shift calculations as listed in Table IXlTI to extract 
6 (r 2 ) 6 ' 7 from the isotope shifts in the Di and D2 lines reported in the literature. The re- 
sults are summarized in Table Kill and also depicted in Fig. [19] in chronological order. The 
large scatter of the data points, much larger than most of the indicated error bars, is the 
first striking impression. Obviously, for many measurements the systematic uncertainties 
were underestimated. In the following, we will shortly discuss the different experiments and 
discuss possible sources of additional or underestimated systematic uncertainties. 

The most reliable early result is probably the measurement on Li + by Riis et al. 
This experiment was performed using saturation spectroscopy on a fast ion beam in collinear 
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TABLE XI: Measured isotope shifts bu^ in the 2s 2 S\/2 — > 3s 2 Si/2 transition as obtained in 
all beamtimes at GSI and at TRIUMF and correction for the second-order Doppler-shift <5^go£>g- 
All values are in kHz. A set of measurements includes typically 15 to 30 resonances at different 
Ti:sapphire laser powers that were used to obtain one ac Stark shift regression line as shown, 
for example, in Fig. [TTJ Results of sets are listed for the on-line beamtimes at TRIUMF. The 
uncertainty of each set is the la confidence limit of the ac Stark extrapolation. If more than one 
set for an isotope was obtained during a beamtime, the average is listed in addition and printed in 
bold. The weighted average of all measurements for one isotope is marked "Average (ToPLiS)". 
Here, the GSI measurement of 9 Li was not included for reasons discussed in the text (Section IIVD|) , 
The measured isotope shift has to be corrected for the differential second-order Doppler shift ^gQ^g 
according to Eq. (|77|) . The final corrected values are marked "Corrected (ToPLiS)". Additional 
systematic uncertainties of 20 kHz + A(Su^q DS ) (uncertainty of the second-order Doppler shift 
coorrection) must be considered for the isotope shift of 7 ' 8 ' 9 Li and of As ys t = 85 kHz as estimated 



for the isotope shift of 11 Li and discussed in detail in Section IIVD1 The total uncertainty of the 
final corrected value is the geometrical sum of the statistical and the systematic uncertainty. 



Facility Year 


bv^ 


^6,8 


Sv 6 > 9 


^6,11 


GSI 2003 


11453 950(130) 


20 089 705(112) 


26 787 087(175) a 




TRIUMF 06/2004 


11453 973(62) 








TRIUMF 09/2004 
Set 1 
Set 2 
Set 3 


11453 982(35) 

11454 058(63) 
11453 993(56) 
11453 889(64) 


20 089 754(81) 


26 787 217(61) 

26 787 363(170) 
26 787 301(105) 
26 787132(82) 




TRIUMF 10/2004 
Set 1 
Set 2 
Set 3 


11453 990(43) 

11453 881(89) 

11454 023(49) 


20 089 767(57) 

20 089 685(104) 
20 089 705(157) 
20 089 826(76) 


26 787 297(55) 


36 555 210(76) 


Average (ToPLiS) 
OU SODS 

Corrected (ToPLiS) 


11453 981(24) 
11(4) 

11453 970(34) 


20 089 754(43) 
19(6) 

20 089 735(50) 


26 787 261(41) 
25(8) 

26 787 236(50) 


36 555 210(66) 
34(11) 

36 555176(108) 


PNNL [25J 


11453 734(30) 









a Not included in the total average for reasons discussed in Section IIVD1 
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FIG. 18: (Color online) Results of isotope shift measurements in the 2s — > 3s transition for 7 Li, 



8 Li and 9 Li at TRIUMF (left) as listed in Table ED The average of all TR 
(A) is plotted on the right and compared to the results of the GSI beamtime 



UMF measurements 
l| f or 7,8,9 Li ^ and 



a measurement of the stable isotopes (v) at the Pacific Northwest National Laboratory (PNNL) 



with a different technique 



251 ] . Uncertainties include the statistical and the systematic uncertainty 



as listed in Table IXj and discussed in the text. 
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and anticollinear geometry to obtain absolute transition frequencies in the ls2s 3 Si — > 
ls2p 3 Po,i,2 transition for the two isotopes. The difference was then calculated to obtain the 
isotope shift. A crossed-beam experiment 89j, where the laser crossed the ion beam under 
90°, was in contradiction for the ls2p 3 Pi_2 fine structure splitting but the result by Riis 
et al. was confirmed in a later experiment using electrooptic modulation 90|. The change 
in charge radius between the stable isotopes extracted from the Li + measurements [l4] is in 
very good agreement with the elastic electron scattering result. 

In contrast, measurements of isotope shifts in the Di and D2 lines show a rather wide 
range of results, which must be attributed to a variety of uncorrected or underestimated 
systematic errors. For example, a recent atomic beam measurement using diode laser res- 
onance fluorescence spectroscopy 91| reports isotope shifts in the Dx and D 2 lines, and 
claims to have determined the absolute transition frequencies with accuracies of 30-60 kHz, 
yet the results for the two lines are internally inconsistent, if the general theories of iso- 
tope shifts and hyperfine structure are to be believed. This claimed precision is much less 
than 1% of the observed linewidth, in a one-photon transition subject to residual Doppler 
shifts and optical pumping distortions. This becomes even more challenging for the D2 
transition, where the hyperfine structure is only partially resolved. Most authors fit the ob- 
served structure by fixing the hyperfine structure interval and intensity ratios according to 
Clebsch-Gordan coefficients. However, optical pumping effects can vary the component line 
intensities dramatically, even at excitation intensities well below saturation 92|, |93| . Further, 
linear polarization and magnetic shielding 9jJ] do not prevent optical pumping; rather one 
should use vanishingly small excitation intensity and/or (pseudo) random polarization. In 
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95| it was tried to account for such effects in the D 2 transition by simulating the effect 



of optical pumping. A shift of almost 2 MHz is proposed for one of the hyperfine structure 
components, but an uncertainty estimation for this rather large correction is not given. 

There are two additional other points which must be considered carefully and may be rel- 
evant for all works on the two lithium resonance lines (Di^): Firstly, exact perpendicularity 
between the atomic beam and the laser is needed to avoid differential Doppler shifts between 
different isotopes. This is particularly important for light-massed isotopes and most works 
make specific efforts to minimize this uncertainty. Nonetheless, at 600 K a misalignment of 
1 mrad from perpendicularity between laser and atomic beam leads to a differential Doppler 
shift of 70 kHz. Secondly, all measurements to date on the D 12 lines uses interferometers for 
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essentially a wavelength measurement. At the level of precision needed for measurements 
of isotope shifts in lithium, interferometric measurements are limited by systematic errors 
associated with pointing errors and mirror phase shifts. As discussed above, resulting scale 
uncertainties are most probably the reason for the deviation of our early 2s — > 3s measure- 
ment and may also contribute to the Di measurement. In comparison, all measurements 
in the 2s — > 3s transition reported here, are based on direct frequency measurements and 
do not suffer from scale errors. There are only two sets of D12 measurements that result 
in consistent (5r 2 ) for both transitions. These are the measurement by Sansonetti and 
coworkers in 1995 96J, with relatively large uncertainties, and a more recent measurement 
by Noble and coworkers in 2006 with considerably smaller uncertainty estimate. The 
latter group previously reported results 95j using the same technique that are in relatively 
poor agreement: The Di isotope shift is shifted by 2a whereas the D2 line is shifted by 4a. 
In Ref. 94j the authors disown their value for the Do and claim that the previous work had 



insufficient resolution. Other measurements [91|, |95|, |97|] yield results that are inconsistent 
both internally, with respect to 5 (r 2 ) values derived for the two transitions, and with each 
other. Hence, also the splitting isotope shift disagrees with theory. 

It is obvious that an accurate determination of the isotope shifts in the D x and D 2 lines of 
lithium is very difficult. The Li + measurements [3] have the clear advantage that a Doppler- 
free technique was used which provides the required accuracy. Similarly, measurements in 
the 2s — > 3s transition, as those reported here, are also to first-order Doppler-free and are 
not hampered by unresolved hyperfine splitting. 

The first isotope shift measurement on the 2s 2 Si/2 —> 3s 2 S\/2 transition reported by 



some of us in 



251 ] was discussed above and likely suffers from a scale error not accounted for in 



the original estimation of systematic uncertainties, because it deviates in the Di transition, 
similarly as in the two-photon transition, from more recent and more accurate measurements. 
The experiments on the different transitions in lithium have also been summarized and 
critically evaluated in [98]. 

In summary, some recent measurements on the Di and D2 lines resulted in a conclusive 
picture for 5 (r 2 ) 6 ' 7 : The optical results reported in 14, 9^] and in this work (which includes 
|l|, Q]) cover all discussed transitions and are consistent within their uncertainties. The 
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FIG. 19: (Color online) Experimental results for the difference in the mean square charge radius 
5 (V 2 ) 6 ' 7 between the stable isotopes 6 Li and 7 Li. Besides the elastic electron scattering (square), 
all other experimental data are obtained from optical isotope shift measurements in different tran- 
sitions in Li + and neutral lithium as indicated in the figure legend and discussed in the text. 
Experimental results are depicted in chronological order. The simple average of the consistent 
data from the individual atomic lithium transitions (solid symbols) is indicated by the solid line 
with the respective standard deviation represented by the dashed lines. All values and respective 
references are listed in Table IXIIl For more discussions see text. 

simple average 5 of 5 (r 2 ) 6 ' 7 obtained from these optical investigations is -0.734(40) fm 2 and 
agrees with the much less accurate result obtained by elastic electron scattering. Thus, 
we conclude that the optical isotope shift measurements combined with atomic mass-shift 



5 Since there are obvious difficulties with the uncertainty estimations at least in some of the previously 
published results, we prefer not to use the uncertainty estimations for a weighted average across all 
transitions. 
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TABLE XII: Changes in the mean square nuclear charge radius 5(r 2 ) 6,7 between the stable isotopes 
6 ' 7 Li from elastic electron scattering and optical isotope shifts in different transitions. Charge radii 



from optical isotope shifts were recalculated using the experimental value Su^ of the respective 
reference and the latest mass shift calculations (D\ t 2 transitions from 5s|, 2s — > 3s from this work). 
Mass shifts <5i^g and field shift coefficients Cq : a are listed for each transition only once. For the 
three transitions in neutral lithium the average values and the corresponding standard deviation is 
listed. The overall average of all measurements is -0.63(20) fm 2 , where the uncertainty is given by 
the standard deviation. The average of the bold printed <5(r 2 ) 6 ' 7 values from [ljlQ] and this work 
is -0.734(40) fm 2 and is indicated in Fig. 1191 This subset was chosen since they give internally and 
mutually consistent values for 5(r 2 ) 6,7 as discussed in the text. 



Transition 


i 6,7 
fe Exp 


c 6,7 


C6,A 


<5<r 2 ) 6 ' 7 


Year & Ref. 




(MHz) 


(MHz) 


(MHz/fm 2 ) 


(fm 2 ) 




Electron scattering 


- 






-0.79(29) 


1972 [84-861 


Li + ls2s 3 Si -> ls2p 3 P 


34 747.73(55) 


34 740.17(3) 


9.705 


-0.78(6) 


1994 [14] 


Li + ls2s 3 Si -> ls2p 3 Pi 


34 747.46(67) 


34 739.87(3) 


9.705 


-0.78(7) 


1994 [14] 


Li + ls2s 3 Si ->■ ls2p 3 P 2 


34 748.91(62) 


34 742.71(3) 


9.705 


-0.64(6) 


1994 [14] 


Weighted Average 








-0.735(36) 


1994 [141 


Li 2s 2 S 1/2 -> 2p 2 P 1/2 


10 532.9(6) 
10533.13(15) 
10 534.26(13) 
10 533.160(68) 
10 534.039(70) 
10 534.215(32) 


10 532.111(6) 


-2.457 


-0.32(24) 
-0.415(61) 
-0.875(53) 
-0.427(28) 
-0.785(29) 
-0.857(16) 


1995 [96] 

1996 [9JJ 
2003 [95] 
2003 [25] 

2006 [94] 

2007 [91] 


Average 


10 533.62(62) 






-0.61(25) 




Li 25 2 S 1/2 -> 2p 2 P 3/2 


10 533.5 (5) 
10 534.90(15) 
10 533.59(14) 
10 534.194(104) 
10 533.352(68) 


10 532.506(6) 


-2.457 


-0.41(20) 
-0.975(61) 
-0.441(57) 
-0.687(42) 
-0.344(28) 


1995 [96] 

1996 [9JJ 
2003 [95] 

2006 [94] 

2007 [91] 


Average 


10 533.91(64) 






-0.57(26) 




Li 2s 2 Si /2 -+3s 2 S 1/2 


11453.734(30) 

11453.95(13) 

11453.970(34) 


11452.821(2) 


-1.572 


-0.581(19) 
-0.719(83) 
-0.731(22) 


2003 [25] 

2004 [1, 2] a 
2010 this work 


Average 


11453.88(11) 






-0.673(71) 





a Result included in the final value of this work. 
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calculations can provide a consistent picture for all transitions investigated so far. However, 
the situation in the D lines is far from being satifactory for both, the isotope shift and the 
splitting isotope shift, due to the large fluctuations in the results. 

4- Isotope shifts of 8,9 Li 

The isotope shifts of the isotopes 8 Li and 9 Li relative to 6 Li were determined similar to 
the measurements of the stable isotopes. A peculiarity of the radioactive isotope scans is 
the fact that ions implanted into the channeltron detector will inevitably decay there and 
produce /3-delayed a particles which will again trigger channeltron signals. This will on 
average result in more than one event per ionized atom and, thus, increase the detection 
efficiency However, one has to ensure that the additional events are counted in the channel 
of that frequency applied to produce the photo ion. Otherwise, the line profile will be skewed 
depending on the scanning direction and the resonance center will be shifted. To avoid this 
effect, the radioactive ion beam was turned on only for 10 s for each frequency channel and 
turned off thereafter for at least five half-lives of the respective isotope (5 s for 8 Li, 1 s 
for 9,11 Li). The events occurring during the beam-off time were still registered and added 
into the corresponding frequency channel. This was possible because the channeltron had 
a very small dark count rate on the order of 10 mHz and the applied detection technique 
of resonance ionization mass spectrometry is, apart from the dark count rate, practically 
background-free. We also looked for events from /3-decay of the n Li daughter n Be (T1/2 = 
13.8 s) and found nothing above background." 

The intensity of the radioactive ion beam was monitored with a Nal detector and a 
plastic scintillator mounted outside the vacuum chamber but close to the carbon catcher 
foil, detecting the 7-rays produced by (3 Bremsstrahlung in the stainless steel chamber and 
/9-delayed neutrons. 

Typical resonance profiles of 8 Li and 9 Li are shown in Fig. [20] (a,c). The count rate in 
each channel iVR aw was normalized to changing Ti:sapphire laser powers and to variations 
in the ion beam intensity according to 

*»»=Mwr-(<^>r <79) 

N y is the 7-ray intensity observed with the plastic scintillator, and Pxiisapphiro is the 
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Ti:sapphire laser power recorded with the photodiode behind the resonator while record- 
ing the respective channel. (iV 7 ) and (-PTi:sa PP hire) are the average values during the com- 
plete scanning period. A correction for the 7-ray background was not necessary for 8 ' 9 Li 
since the count rate with the ion beam turned on was several orders of magnitude larger 
than the background with the ion beam turned off. The error bars of the data points in 
Fig. [20] include the statistical uncertainty as well as the propagated uncertainty from the 
normalization procedure. 

Fitting was performed as for the stable isotopes and the e.g. was calculated from the 
hyperfme component positions according to Eq. ( 16T)|) . As in the case of the stable isotopes, 
a small asymmetry in the line profile is observable in the residuals. It appears that the 
asymmetry is inverted for 9 Li compared to the stable isotopes and 8 Li as is the case for 
the order of the hyperfme components. Both is caused by the inverted frequency scale for 
the heavier isotope 9 Li. The frequency of the Ti:sapphire laser is measured relative to the 
iodine-locked diode laser (z^dl)- Resonance frequencies for 6,7 ' 8 Li are below the diode laser 
frequency while the 9 ' n Li resonances appear at higher frequencies. 

Results of ac Stark shift measurements for 8 ' 9 Li are shown in Fig. EH] (b,d). The seemingly 
opposite slope of the lines is again due to the inverted frequency scale. As for the stable 
isotopes, the isotope shift is determined from the intercept of the regression line with the 
y-axis referenced to a similar measurement of the reference isotope 6 Li. Measurements of 
6 Li were regulary interspersed between the ones on short-lived isotopes. The results for 
the isotope shift of 8 Li and 9 Li are included in Table IXII Four sets of measurements were 
obtained for each isotope during the two beamtimes at TRIUMF. The results of the different 
measurements and their averages are plotted in Fig. [18] Systematic uncertainties are similar 
as for the stable isotopes. The comparison with the measurements from GSI (♦) shows good 
agreement for both isotopes. 

5. Isotope Shift of n Li 

The very short-lived isotope n Li was available at TRIUMF with a yield of about 35 000 
ions/s. A single scan of a n Li resonance curve was acquired within approximately 10 min. 
Typically 5 to 7 consecutive individual scans were added to obtain a resonance signal as 
shown in Fig. [21] The small production rate of n Li as compared with the other isotopes 
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FIG. 20: (Color online) Typical spectra (a,c) and ac Stark shift results (b,d) for 8 Li and 9 Li. Data 
points of the spectra are normalized to intensity variations of the ion beam and laser beam according 
to Eq. (|79p . Voigt profiles with Doppler pedestal are fitted to the hyperfine resonances with the 
same procedure as used for the stable isotopes. The solid lines in (a) and (c) are the common fit 
results while the dashed lines are the line profiles of the individuals components. Results of ac 
Stark shift measurements are plotted in (b) and (d) for 8 Li and 9 Li, respectively. The solid line 
is the result of a linear regression and the dashed lines indicate the 95% confidence band. The 
intercept of the regression lines with the y-axis is used to determine the isotope shift. 
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did not allow for full Thsapphire laser power-series measurements since it was not possible 
to collect sufficient statistics at low laser power within a reasonable time. Thus, a different 
approach was used to correct for the ac Stark shift of this isotope: All n Li spectra were 
recorded at maximum laser power of the Ti: sapphire laser and the average power value during 
the scan was evaluated. The count rate was normalized for ion beam and laser intensity 
fluctuations as described above for 8 ' 9 Li, but this time the background in the 7 and neutron 
count rate of the plastic scintillator had to be taken into account. Therefore, the 1 s pause 
after each frequency step was used to determine the ion-beam independent background and 
to subtract it before normalization. The resonance signal was then fitted with the same line 
profile as applied to the other isotopes and the hyperfine e.g. was determined. For reference, 
a complete ac Stark shift power-series measurement of 6 Li was performed at least once per 
day. Then, the ac Stark shift function determined for 6 Li 

K 6 Li) = z/ ( 6 Li) + 6J*-ppU» . P Ti:sapphire (80) 

was evaluated at the average power (-PrusapphmJuLi during the n Li measurements and the 
isotope shift was determined according to 

Su 6 ' 11 = u( n U) - (^ ( 6 Li) + fcffc-w**- . (P Ti;sapphirc )n Li ) . (81) 

This procedure requires that the ac Stark shift of both isotopes have the same slope 
frTi.sapphire. This is theoretically expected since all resonance frequencies in the 2s — > 3s 
transition are far-off from any allowed El transition that would connect the 2s or 3s level 
to a p state. In fact, this was also experimentally verified for all other lithium isotopes as 
will be described in the next Section. 

In total, 24 resonances could be measured for 11 Li. The distribution of the extracted 
isotope shift values Si; 6,11 is shown in Fig. [22] (open circles). The weighted average is included 
in Table IXII The uncertainty was estimated as the standard error of the mean as depicted 
in the figure. This is in accordance with the uncertainties obtained for the more abundant 
isotopes discussed above. 

D. Systematic uncertainties in the isotope shift 

Several effects that can possibly affect the extracted resonance position were investigated. 
These include deviations from the Lorentzian, respectively Voigtian resonance lineshape, 
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FIG. 21: (Color online) Typical spectrum of n Li as it was obtained in the beamtime at TRIUMF. 
Data points are normalized to intensity fluctuations of the ion and laser beams. The solid line 
is again the common fit function of the two hyperfine components while the dotted lines are the 
contributions of the individual hyperfine components. 

isotope-dependent slopes or nonlinearities in the ac Stark shift correction, and the influence 
of stray magnetic fields. 

It was found that the most crucial effect is an isotope-dependent slope of the regression 
line in the ac Stark correction. This does not matter for the stable and the less exotic 
isotopes 6_9 Li since those were all measured at different powers and extrapolated to zero 
laser intensity. However, the slope for 11 Li could not be measured experimentally and a 
difference in the slope between 11 Li and the reference isotope 6 Li would directly influence 
the extracted isotope shift. Therefore, the slopes from all 6_9 Li ac Stark shift measurements 
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FIG. 22: (Color online) Distribution and the average of the n Li - 6 Li isotope shift (e.g.) of all 24 
n Li resonances obtained at TRIUMF. The dotted lines represent the standard deviation (a) of the 
individual measurements and the dashed ones the standard error of the mean. The average value 
is plotted as solid circle. A histogram of the measured frequencies is depicted on the right. 

at TRIUMF were compared after the first on-line beamtime. Previous measurements at 
GSI showed no indication for an isotope-dependent slope in the ac Stark shift of 6 ' 7 ' 8 Li 
within the level of accuracy of approximately 2 kHz/mW. However, at TRIUMF the ac 
Stark shift showed larger isotope-dependent trends. This can be seen in the left part of 
Fig. [23] where this is indicated by the straight lines. There is a clear offset between 7 Li and 
8 Li, and the reference isotope 6 Li shows even a small drift. This effect was found to be an 
artefact caused by interference effects in the viewport behind the optical resonator. The 
geometry of the vacuum chamber, the resonator and the viewport is depicted in Fig. [TOl 
While the entrance window was antireflection coated, the exit window behind the resonator 
was uncoated. The plane-parallel surfaces of the 5 mm thick window gave rise to an etalon 
effect, and thus the amount of light detected on the photodiodes was wavelength-dependent 
and therefore different from isotope to isotope since the free-spectral-range of the window- 
etalon was comparable to the isotope shift between 6 Li and 11 Li. A direct measurement of 
this effect was performed by comparing the amount of light entering the vacuum chamber 
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FIG. 23: (Color online) Experimentally determined ac Stark shift coefficients b 1 1,sapp irc for all 
measurement sets and all isotopes conducted at TRIUMF. The periods of the two beamtimes 
are indicated at the bottom. Measurements 1-37 were performed with a plane parallel viewport 
window while later measurements (gray background) were carried out with a wedged window to 
avoid interference effects. The lines are included to guide the eye along the measurements of the 
respective isotope (left) and to indicate the common drift after the replacement of the viewport 
with a wedged window. The highlighted part on the right indicates the period in which all n Li 
resonances were recorded. 

and leaving it. The etalon effect was clearly observed and the plane-parallel window was 
replaced by a wedged window. A second measurement after the replacement did not show 
a detectable wavelength-dependent variation anymore. 

Ac Stark shift coefficients that were measured after the replacement are plotted in the 
right part of Fig. [23j There is still a slow drift of the coefficients with time, which might 
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be caused by a slow degradation of the mirror reflectivity, but all isotopes show equal co- 
efficients within the uncertainty. To estimate the remaining systematic uncertainty for the 
11 Li measurements, the standard deviation of the slope coefficients for the reference iso- 
tope within the periode of n Li measurements, as indicated in the figure, was determined 
(a = 0.81 kHz/mW) and multiplied with the average Ti:sapphire laser power on the pho- 
todiode during all measurements (30 mW). Since much less measurements were performed 
for n Li than for the other isotopes, the ac Stark shift of the dye laser might not be fully in- 
cluded in the statistical uncertainty. Therefore, a systematic contribution was estimated by 
multiplying the slope of the ac Stark shift caused by the dye laser (Eq. (|71|) . obtained from 
Fig. [To*]) with the standard deviation of the dye laser power as measured on the photodiode 
during all n Li measurements. This resulted in systematic uncertainties of 46 kHz for the 
shift induced by the intensity of the Ti:sapphire laser light and of 8 kHz by that of the dye 
laser. These values are listed in Table IXIIII 

The observation of this interference effect is the reason why the 9 Li isotope shift mea- 
surement at GSI was excluded in the final average: Insufficient statistics did not allow us 
to perform an ac Stark shift measurement for this isotope at GSI. Instead, a measurement 
as for 11 Li at TRIUMF was carried out but with the plane-parallel viewport. Therefore the 
value might be systematically shifted. However, the reasonable agreement with the TRI- 
UMF data within the uncertainty shows that the size of a possible systematic contribution 
was reasonably estimated in [lj]. 

Recoil effects are not important for two-photon transitions, since the momenta of the two 
photons cancel and the velocity of the atom stays constant during the excitation process. 
The differential second-order Doppler shift correction £^sod was already discussed and its 
uncertainty is included in the systematic uncertainty. 

The influence of unresolved Zeeman splitting in the hyperfine structure has been theoret- 
ically estimated by calculating the effect of the 0.7 G stray magnetic field measured at the 
beamline. A large part of the induced shifts in the 2s and 3s levels cancels. The remaining 
part shifts the peak positions of the two hyperfine components by only about 10 kHz in 
opposite directions. Thus, the e.g. of the hyperfine structure is expected to stay constant. 
For safety, the total expected shift of ps 10 kHz is included in the error budget for the isotope 
shift (Table IXIj) as well as for the hyperfine splitting (Table |X|) . 

Finally, the small asymmetry in the peak profile has to be considered. This was observed 
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27] , and from second-order 



TABLE XIII: Summary of statistical and systematic uncertainties of the corrected n Li - 6 Li isotope 
shift 6i/jg . Systematic uncertainties arise from slope variations in the ac Stark shift induced by 
the Ti:sapphire or dye laser induced between n Li and 6 Li, unresolved Zeeman splitting caused 
by the stray magnetic fields (0.7 G), a lineshape asymmetry caused by a position-dependent ac 
Stark shift in the laser beam intensity profile as discussed in detail in 
Doppler shift. All systematic contributions are linearly added to A syst . The total uncertainty is 
the geometrical sum of the statistical and the systematic contributions. 

Effect Contribution 

Statistical uncertainty 

Standard error of the mean for 24 measurements 66 kHz 

Systematic uncertainty 

ac Stark shift induced by the Ti:sapphire laser 46 kHz 

ac Stark shift induced by the dye laser 8 kHz 

Unresolved Zeeman splitting 10 kHz 

Lineshape asymmetry 10 kHz 

Second-order Doppler shift 11 kHz 

Total A Syst 85 kHz 

RMS Total 109 kHz 

much more pronounced in off-line measurements at GSI which were carried out very recently 
in order to determine the absolute frequency of the 2s — > 3s transition. The asymmetry was 
attributed to the position-dependent ac Stark shift that the atoms experience when they 
cross the Gaussian laser beam profile . The strongest shift occurs in the beam center but 
all atoms that are excited and ionized in the wings of the laser beam are less ac Stark shifted 



and are therefore causing an asymmetry of the lineshape. 



his assumption was supported 



3- 



There it was shown that 



by lineshape simulations which are presented in detail in 
the asymmetry leads to a small shift in the position and results in a small deviation of the 
ac Stark shift behavior from a straight line at very low power. This is crucial concerning 
the absolute transition frequency, but the isotope dependence was found to be very small 
(< 10 kHz). This value was added to the systematic uncertainty for all isotopes. 
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All effects contributing to the uncertainty of the n Li - 6 Li isotope shift are summarized 
in Table IXIII1 Even though the systematic uncertainties are mutually independent and will 
not necessarily cause a shift in the same direction, we have added them linearly to obtain the 
total systematic uncertainty which is then geometrically added to the statistical uncertainty 
In total, a relative uncertainty of the 6 ' n Li isotope shift of 3 • 10~ 6 was obtained in this work. 

For the other isotopes, power fluctuations of the dye laser are included in the statistical 
uncertainty given in Table IXll The Zeeman splitting (10 kHz), the lineshape asymmetry 
(10 kHz) and the uncertainty in the second-order Doppler shift are correlated for all isotopes 
and have therefore been linearly added to the total systematic uncertainty. The geometric 
sum of the systematic and the uncorrelated statistical contribution is then taken as the final 
isotope shift as listed in Table IXII 



V. CHANGES IN THE MEAN SQUARE CHARGE RADIUS 



The measured isotope shifts (Table IXip can now be combined with the mass shift cal- 
culations discussed in Section III Fl (Table IHIj) to extract by use of Eq. (|62|) the change in 
the rms nuclear charge radii 5 (r 2 ) along the lithium chain of isotopes. Results are listed in 
Table IXIVI and plotted in Fig. [2H relative to 6 Li. This table represents the summary of the 
complete analysis of the experimental results and theoretical mass shift calculations. We 
rat slightly varying mass shift values along previous publications on this 
are not caused by a model-dependence of the mass shift calculations 



want to stress t 
topic 1 Q, 



57, 



OS 



but several improvements: First, since the start of the experiments the accuracy of the 
atomic calculations could be improved by an order of magnitude; secon d, m ore accurate 



atomic masses became available by mass measurements at MISTRAL 



99 



100| and TITAN 



53[| ; third, the influence of the nuclear polarizability on the transition frequency could be 



calculated for the first time 



29] . Hence, previous values should not be used anymore for 



calculating the nuclear charge radii. 

Figure [24] clearly depicts the trend in the nuclear charge radii: The size of the proton 
distribution shrinks continuously from 6 Li to 9 Li. As already stated in js| this reduction is 
attributed to the strong clusterization of the lithium nuclei. Clusterization is a phenomenon 
that often appears in light nuclei; they tend to behave as composed of well-defined subsys- 
tems, called "clusters." Experimentally, this is reflected in the fact that a-particles, tritons 
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FIG. 24: (Color online) Changes 5(r 2 )®' A in the mean square charge radius along the lithium isotope 
chain with respect to the mean square nuclear charge radius of 6 Li. The values are obtained from 
the measured isotope shift (Table IXI|) and theoretical mass shift calculations (Table IIII|) using 
Eq. & 



(t = 3 H), and helions ( 3 He) can be knocked out from a number of nuclei with substantial 
probability. Often a clusterization is observed even in heavier stable nuclei along the N = Z 
line, for example in 12 C or 16 0. According to nuclear structure calculations, 6 Li can be de- 
scribed as a + d, 7 Li as a + t, and the heavier isotopes as a + t + xn with x = 1, 2 for 8 ' 9 Li, 
respectively. The consequence of clusterization is a strong center-of-mass (cm) motion of 
the clusters that leads to larger proton distributions. This cm motion is obviously reduced 
with increasing number of neutrons, hence, the nuclear charge radius decreases. This trend 
is reversed when going from 9 Li to n Li. Here the charge radius increases again and n Li is 
similar in size to 7 Li. There are in general two effects that can be made responsible for this 
increase: Either the 9 Li core of n Li is changed in its internal structure, i.e. core polariza- 
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TABLE XIV: Field shifts contributions Su F ' s in the 2s — > 3s transition of the lithium isotopes 
(extracted from the measured isotope shifts in Table IXll by subtracting the theoretically evaluated 
mass shift contributions listed in Table IIIip , field shift coefficients C as defined by Eq. ()60p , and 
the change in the mean square nuclear charge radius 5(r 2 ) e ' A between the isotope ^Li and the 
reference isotope 6 Li calculated using Eq. (|62[) . 



A Li Field Shift (MHz) C (MHz/fm 2 ) 5(r 2 ) 6 > A (fm 2 ) 

7 Li 1.149(34) -1.5719(16) -0.731(22) 

8 Li 1.933(51) -1.5719(16) -1.230(32) 

9 Li 2.615(50) -1.5720(16) -1.663(32) 

n Li 0.852(108) -1.5703(16) -0.543(69) 



tion takes place, or the correlation of the halo neutrons induces a cm motion of the 9 Li core 
nucleus. 

Core polarization is related to the neutron- 9 Li interaction and the internal structure of 
9 Li, whereas the center-of-mass motion gives insight into the neutron-neutron correlation in 
the three-body system. Pure three-body calculations assume the core to be inert ("frozen 
core" approximation) whereas more elaborated nuclear models try to take both effects into 
account. There is experimental evidence that the frozen core assumption is a very good ap- 
proximation, however, both effects might contribute to the total change in the charge radius 
between 9,11 Li and it is impossible to disentangle their contributions without consideration 
urther experimental data or theory. More complete theoretical models were considered 



of 
in 
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3[ and will be discussed in more detail in a following publication 
Comparison with nuclear models can best be performed if absolute charge radii are calcu- 
lated from the measured isotope shift. However, we need the charge radius of one reference 
isotope to fix the absolute charge radius. The value for the charge radius obtained by this 
procedure is model-dependent if the reference radius is obtained using a nuclear model. 
This is the case for the charge radii from elastic electron scattering of lithium, because the 



model-independent approach gave rather large uncertainties 



86]. Since the value of 5 (r 2 ) 6 ' 7 



obtained in our isotope shift meaurement is slightly different from that obtained by elastic 



electron scattering 



87j (even though within uncertainties), the absolute charge radii differ 



if they are based either on the charge radius reported for 6 Li 



85 



86[ or that reported for 
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84j . In previous publications, we have used r c ( 7 Li) = 2.39(3) fm from 84( which is 



not model-independent and a recent analysis of the world scattering data revealed that the 
charge radius of 6 Li is better understood. Since this topic is of major importance for the 
comparison with nuclear structure models and for conclusions about the structure of the 



halo nucleus, it is discussed in more detail in the following publication i 



VI. SUMMARY 



This long writeup reports about the isotope shift measurement of the complete lithium 
isotope chain from the stable isotopes up to the exotic two-neutron halo nucleus 11 Li. The 
theory of atomic mass shift calculations is explained in detail and the most accurate results 
for the mass shifts in the 2s 2 Si/2 —> 3s 2 Si/2 two-photon transition are presented. The 
experimental setup for high-accuracy resonance ionization mass spectroscopy with cw lasers, 
including a two-photon transition is described with all important details to judge the quality 
and accuracy of the data obtained. All sources of possible systematic uncertainties in these 
measurements that are known to us are evaluated and discussed. Compared to previous 
publications, a small additional correction for the second-order Doppler shift is included 
in the isotope shift result. Combination of the accurate measurements with the mass shift 
calculations allows us to extract the changes in the nuclear charge radii along the isotope 
chain spanning from the stable 6 Li to the halo nucleus n Li. Absolute nuclear charge radii will 
be extracted and compared with predictions by nuclear models in a following publication. 
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